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FOREWORD 


This  report  represents  one  phase  of  an  analytical  investigation  of  handling 
qualities  for  multiloop  airframe-hunan  pilot  systems.  The  research  reported 
was  sponsored  by  the  Flight  Control  Laboratory,  Aeronautical  Systems  Division, 
Air  Force  Systems  Command,  as  part  of  Project  No.  8219,  Task  No.  821 905.  It 
was  conducted  at  Systems  Technology,  Inc.,  under  Contract  No.  AF  35(6l6)-8024. 
The  ASD  project  engineer  was  Mr.  R.  J.  Wasicko. 


The  authors  would  like  to  acknowledge  the  fine  work  of  J.  Taira  and 
Dolores  Lewis  in  the  preparation  of  the  final  manuscript. 
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ABSTRACT 


Two  related  facets  of  feedback  system  analysis  are  considered — the 
calculation  of  closed- loop  response  and  the  determination  of  the  effects 
on  closed- loop  behavior  of  variations  in  open-loop  parameters .  The  con¬ 
nection  between  the  modal  response  coefficients  (partial  fraction  expan¬ 
sion  coefficients  of  the  closed- loop  transfer  function) ,  the  sensitivity 
of  the  closed- loop  poles  to  variations  in  the  open-loop  gains,  poles,  and 
zeros,  and  classical  sensitivity  are  developed.  A  con^jrehensive  summary 
is  given  of  methods  for  deteimining  modal  response  coefficients  and 
sensitivities  from  open-  or  closed-loop  transfer  function  representation. 
Response  formulas  for  periodic  or  power  series  inputs  are  derived  in 
terms  of  the  modal  response  coefficients,  with  conventiorial  error  coeffi¬ 
cients  as  a  special  case.  The  initial  developments  are  for  single-loop 
systems  with  first-order  closed- loop  poles;  these  are  extended  to  cover 
multiple-order  closed- loop  poles,  and  generalized  to  multiloop  systems. 
Examples  are  given  using  jto,  I,  cr,  and  shifted  Bode  plots,  root  locus, 
and  analytical  transfer  function  representations. 
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SYMBOIS 


Real  part  of  z  or  p 

Coefficient  of  in  numerator  of  G(s)  when  G(s)  is  in  root-locus  form 
Amplitude  of  cos  o^t  input 
Amplitude  of  G(s) 

Bode-form  numerator  of  G(s) 
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Coefficient  of  in  denominator  of  G(s)  when  G(s)  is  in  root-locus 

form 

Amplitude  of  sin  o^t  input 
Bode-form  denominator  of  G(s) 

System  output 

ySITbg 

Distance  imaginary  axis  is  shifted  to  the  left 
Coefficient  of  t^/MI  input 
Error,  r(t)  -  c(t) 

Error  coefficient 

Numerator  of  closed- loop  transfer  function 

Denominator  term  of  closed- loop  transfer  function  (denominator  is  1 +  F2) 
Open-loop  transfer  function 
Controller  transfer  function 
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Feedback  transfer  function 

Transfer  function  of  kth  element  in  multiloop  system 
Closed- loop  output/input  transfer  function 
Closed- loop  error/input  transfer  function 
Altitude 
Index 

Altitude/elevator  transfer  function 
Index 

System  response  to  unit  step  input 
Index 


Index 
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Bode  gain  of  G(s) 

Number  of  transfer  functions  in  multiloop  system 
Inverse  Laplace  transform 

Number  of  poles  of  G(s)  minus  number  of  zeros 
Index  for  power  series  inputs 
Number  of  zeros  of  G(s) 

Order  of  closed- loop  pole 


Noise 


Q*(s) 


Negative  of  open- loop  pole 
Negative  of  closed-loop  pole 
Modal- response  coefficient 
See  Eq  8,  9^^  1 66 
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SECTION  I 


INTRODUCTION 

A.  GENERAL  LINEAR  SYSTEM  ANALYSIS  PROBLEM 
1 .  The  Analysis  Problem  Defined 


In  this  report  detailed  consideration  is  given  to  two  major  facets  of 
feedback  system  analysis — the  effects  of  component  variations  on  system  transfer 
functions,  and  the  connection  between  tmnsfer  function  characteristics  and 
system  time  responses.  For  constant- coefficient  linear  systems  typified  by  the 
single-loop  system  shown  in  Fig.  1 ,  the  system  analysis  problem  shall  be  con¬ 
sidered  here  to  consist  of  five  essential  steps: 

a.  Delineation  of  nominal  open- loop  system  characteristics. 

This  step  ordinarily  starts  with  differential  equations 
describing  the  nominal  controlled-element  characteristics 
and  one  or  more  controller  possibilities.  The  step  is 
concluded  when  one  or  more  nominal  open-loop  transfer 
functions,  G(&),  are  available,  in  factored  form,  for 
further  analysis. 

b.  Determination  of  nominal  closed-loop  transfer  functions, 

Gj,(,(s)  and  Gj,g(s),  from  the  open- loop  transfer  function(s). 

c.  Calculation  of  closed-loop  system  time  responses  for 
pertinent  inputs. 

d.  Determination  of  the  changes  in  G(s)  resiilting  from  the 
expected  variations  in  the  controller  and  controlled- 
element  characteristics . 

e.  Consideration  of  the  effects  of  open- loop  system  variations 
on  closed- loop  behavior. 

The  major  topics  covered  in  this  report  are  concerned  with  steps  c  and  e, 
although  desirable  techniques  to  accomplish  steps  c  and  e  are  somewhat  dependent 
on  the  methods  used  in  steps  a  and  b. 


Original  manuscript  submitted  August  1962;  revised  manuscript  released  by 
authors  November  1962  for  publication  as  an  ASD  Technical  Docvunentary  Report. 
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Figure  1 .  Single-Loop  Linear  Feedback  System  and  Basic  Notation 


Delineating  the  open-loop  characteristics  in  terms  of  a  transfer  function 
G(s)  is  extremely  simple  for  time-invariant  linear  systems  because  transform 
methods  can  he  used  to  convert  the  system  differential  equations  to  algebraic 
equations.  This  conversion  permits  the  intermediate  steps  in  an  analysis 
sequence  (e.g.,  reduction  of  simultaneous  equations  using  Cramer’s  rule, 
transfer  fimction  development  and  manipulation)  to  be  carried  out  using 
algebraic  forms.  Most  such  forms  are  rational  polynomials,  and  the  remainder 
can  usually  be  approximated  by  this  type  of  function.  Thus,  with  the  possible 
exception  of  polynomial  factoring  (which  also  enters  step  b),  the  delineation 
of  open- loop  transfer  functions  in  the  form  indicated  in  Fig.  1  is  basically 
elementary,  as  is  the  extension  to  similar  open- loop  transfer  functions  in 
multiloop  systems. 

The  second  step  in  the  analysis  sequence,  i.e.,  given  G(s),  find  Gj.(j(s) 
and/or  Gj.e(s),  is  the  central  problem  of  feedback  systems  analysis.  Trivial  as 
it  may  seem,  amounting  only  to  finding  roots  of  1  +  G(s)  =0  when  G(s)  is  given, 
a  great  deal  of  effort  has  been  devoted  to  finding  better  methods  for  performing 
this  operation.  An  eclectic  view  of  this  step,  for  single-loop  systems,  is 
presented  in  Ref.  1 .  Various  representations  of  the  open-loop  transfer  function 
(such  as  G(ja)),  G(±a),  and  G(|,p)  Bode  plots  and  pole-zero  locations)  are  used 
there  as  part  of  a  unified  combination  of  methods  to  find  closed- loop  transfer 
functions  in  factored  form.  The  techniques  of  Ref.  1  have  also  been  extended  to 
multiloop  systems  in  Ref.  2, 

2.  System  Time  Response  Behavior 


When  closed-loop  transfer  functions  are  available  in  factored  form  step  c 
is,  in  principle,  an  elementary  inversion  process  from  the  transform  to  the  time 
domain.  The  major  operation  involved  can  again  be  algebraic,  e.g.,  the  system 
response  to  a  unit  impulse  input  requires  the  resolution  of  the  transfer  function 
Gj.c(s)  into  a  partial  fraction  series.  This  algebraic  operation  is  then  followed 
with  term-by-term  inversion  to  the  time  domain.  Thus,  using  the  Gj,(>(s)  form 
shown  in  Fig.  1  (for  m>  l). 
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The  first  step  is  to  resolve  Gj,(j(s)  into  partial  fractions  which,  assuming  that 
there  are  no  repeated  closed-loop  poles,  results  in 


m+n  Q. 

Grc(s)  =  i  TT~c 


i=l 


Then,  inverse  transforming. 


-1  m+n  . 

£  prc(-)]  =  E  V 


gives  the  system  time  response  to  the  unit  impulse  ii5)ut. 

Standard  procedures  for  finding  the  modal  response  coefficients — the  Qj^’s 
above — are  routine  but  tedious.  Techniques  for  finding  the  modal  response 
coefficients  from  the  open-loop  transfer  function  6(s)  are  available  but  not 
too  well  known.  So,  the  analyst  often  stops  when  Gj.(.(s)  is  known  in  factored 
form,  secure  in  the  knowledge  that  the  time  variation  portions  (the 
terms)  of  the  system  response  are  defined  by  the  qj|^’s  appearing  explicitly  in 
the  closed- loop  transfer  functions,  and  that  the  relative  magnitudes  of  these 
time  variations  (the  Q^'s)  are  indicated  there  implicitly.  Much  of  this  report 
is  devoted  to  developments  and  methods  which  provide  the  modal  response  coeffi¬ 
cients  directly  from  the  open-  and/or  closed-loop  transfer  function  representa¬ 
tions  used  in  the  second  step  of  the  general  analysis  procedure.  Using  these 
methods,  the  process  of  finding  modal  response  coefficients  (and  system  time 
responses)  becomes  a  mere  adjunct  to  the  second  step  rather  than  a  separate, 
often  neglected,  stage  in  analysis.  Further,  by  extending  the  techniques  to 
multiloop  systems,  the  third  step  in  the  general  analysis  problem  gains  the 
same  status  of  completeness  as  the  first  two  already  possess. 

3.  System  Sensitivity  to  Parameter  Variations 


One  reason  for  the  fifth  step  (effects  of  component  variations  on  system 
characteristics)  is  obvious — the  system  assumed  in  an  analytical  study  will 
never  precisely  match  the  actual  physical  system,  so  it  is  important  to  know 


h 


the  effects  of  possible  variations.  A  second,  core  subtle  reason  is  that  a 
general  knowledge  of  the  effects  of  paraiceter  changes  can  be  used  as  a  guide  to 
system  codifications  which  would  improve  the  over-all  performance.  For  both  of 
these  reasons  measures  of  closed-loop  system  sensitivity  to  open- loop  system 
parameter  variations  are  an  integral  part  of  the  analysis  problem. 

Kotions  about  system  sensitivity  were  in  the  forefront  when  the  feedback 
concept  was  initially  developed.  This  was  natural,  even  uiKivoidable,  since 
feedback  systems  possess  the  "fundamental  physical  pi'operty  that  the  effects  of 
variation  in  the  forward  loop,  whether  they  are  taken  as  changes  in  G{s)  or  as 
departures  from  strict  linearity  or  from  freedom  from  extraneous  noise,  are 
reduced  by  the  factor  1  +  G  in  comparison  with  the  effects  which  would  be 
observed  in  a  non-feedback  system. Accordingly,  sensitivity  measures  were 
indispensable  to  any  rational  discussion  of  feedback  systems,  and  a  useful, 
classical  definition  of  sensitivity  was  made  one  of  the  two  mathematical  defini¬ 
tions  of  feedback  (Ref.  3)-  Except  for  a  minor  modification  (Ref.  4)  this 
definition,  by  Bode  (Ref.  3)^  remained  unchanged  for  over  a  decade.  Perhaps 
this  static  nature,  surrounded  by  dynamic  growth  in  most  other  areas  of  feedback 
systems  engineering,  made  the  concept  fade — in  any  event,  whole  cadres  of 
neophyte  systems  engineers  were  trained  with  precious  little  exposure  to  sensi¬ 
tivity  concepts.  In  recent  years  classic  sensitivity  has  become  a  more  popular 
subject  in  automatic  control  (Ref.  4  and  5) •  Finally,  the  en^hasis  on  pole-zero 
specifications  for  system  characteristics  gave  rise  to  new  conceptions  of  sensi¬ 
tivity,  with  associated  new  measures.  These,  called  here  "gain,”  "(open- loop) 
pole,"  and  "(open-loop)  zero"  sensitivities,  were  evolved  to  account  for  changes 
in  the  position  of  closed- loop  poles  due  to  shifts  or  changes  in  open- loop  gain, 
poles,  or  zeros.  A  number  of  thesis  studies  and  research  efforts,  such  as  those 
reported  in  Ref.  6-13,  have  gone  a  long  way  toward  bringing  this  subject  to  a 
logical  conclusion. 

In  the  simplest  terms,  sensitivity  relationships  connect  open- loop  differen¬ 
tial  variations  with  closed-loop  pole  differential  shifts.  Thus,  the  differential 
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Paraphrased  quotation  from  Ref.  3^  P*  44. 


shift  of  a  closed-loop  pole  is  given  hy 


where  the  factors  are  first-order  sensitivities  and  the  differentials  dzj, 
dn^,  and  dfc  may  he  interpreted  as  incremental  changes  in  the  open-loop  parame- 
ters.*  IThe  ^in  sensitivity,  is  basic,  since  the  pole  and  zero  sensitivities 
can  be  related  to  it,  and  since  it  is  a  coiroonent  of  the  classical  sensitivity. 

When  the  sensitivity  factors  are  known,  they  may  be  combined  with  the  results 
of  the  fourth  step  in  the  general  analysis  sequence  (estimates  for  k,  zj,  and  p^ 
uncertainties)  to  provide  tolerances  on  the  closed- loop  roots.  Further,  a 
knowledge  cf  the  values,  and  of  their  connection  with  system  parameters,  can 
lead  to  system  changes  which  minimize  the  system  variations  to  element  deviations, 

Sensitivities  can  also  be  very  useful  in  the  synthesis  problem.  Once  a 
trial  controller  has  been  formed,  the  sensitivities  permit  the  designer  to 
estimate  what  changes  he  should  make  in  the  controller  to  obtain  a  more  desirable 
set  of  nominal  closed- loop  poles.  This  approach  is  particularly  beneficial  for 
multiloop  systems  where  the  effects  of  inner- loop  parameters  on  the  fint,!  closed- 
loop  poles  are  frequently  obscure. 

Sensitivities  also  have  great  value  when  a  digital  computer  is  used  as  the 
prime  means  of  control  analysis.  By  extending  the  program  to  compute  sensitivi¬ 
ties,  the  designer  can  use  a  manual  or  automatic  iteration  scheme  to  "home  in" 
on  the  best  controller  rather  than  doing  a  costly  parametric  investigation  of 
all  the  controller  parameters. 

On  the  surface,  steps  c  and  e  of  the  analysis  problem  appear  to  be  separate 
and  distinct.  However,  there  is  a  strong  tie  between  the  two — the  gain 


‘Note  that  in  a  strict  mathematical  sense  the  closed- loop  poles  are  -qj, ,  the 
open- loop  zeros  are  -z^,  and  the  open-loop  poles  are  -pj.  The  notation  used  in 
this  report  (Fig.  l)  has  been  chosen  to  agree  with  a  standard  servoanalysis 
convention  which  has  transfer  function  factors  as  (s  +  Zj)  and  (s  +  Pj)* 


sensitivities  being,  in  fact,  equal  to  the  modal  response  coefficients  for 
single-loop  systems  with  single-order  poles.  This  close  connection  provid.es 
ad.ditional  methods  for  the  calculation  of  modal  response  coefficients  and  gain 
sensitivities  from  the  conventional  open-  and  closed-loop  system  representations. 
These  methods  can  also  be  extended  to  multiloop  situations  where  sensitivity 
concepts  should  ultimately  be  very  helpful  in  reducing  to  well-understood  terms 
the  many  puzzling  interactions  present. 

B.  OUTLIlffi  OF  THE  REPORT 

The  remainder  of  this  report  is  organized  into  four  sections  (II  through  V) . 
The  basic  relationships  for  modal  response  coefficients  and  sensitivities  for 
single-loop  systems  are  developed  in  Section  II.  Several  methods  for  evaluating 
the  modal  response  coefficients  and  sensitivities  are  derived  in  Section  III. 
Section  IV  extends  the  developments  to  multiloop  systems,  and  Section  V  summa¬ 
rizes  the  important  results. 

The  developments  of  Section  II  begin  with  the  modal  response  coefficients 
for  single-order  closed-loop  poles.  The  connections  between  the  modal  response 
coefficients  and  the  closed-loop  responses  for  siii5)le  inputs  are  developed. 

Then  the  sensitivities  and  their  relationship  to  the  modal  response  coefficients 
and  classical  sensitivity  are  de..ived.  The  final  article  contains  the  modifica¬ 
tions  and  extensions  to  the  above  which  are  necessary  for  systems  with  multiple- 
order  closed-loop  poles. 

The  evaluation  methods  presented  in  Section  III  are  grouped  as  requiring 

1  .  Direct  calculation 
2.  Root  locus  plots 
5.  Open-loop  Bode  or  |  plots 
4.  Closed-loop  Bode  asymptotes 

The  numerical  results  obtained  in  an  illustrative  example  are  presented  to  indi¬ 
cate  the  accuracies  obtainable  with  each  technique.  Details  of  the  computations 
are  given  in  the  Appendix. 


The  extensions  for  multiloop  systems,  which  are  given  in  Se::tion  IV,  generally 
j.NV-y ■  urallel  the  developments  in  Section  II.  Section  IV  concludes  with  an  example 
for  a  flight  control  system  with  pitch  angle  and  altitude  feedbacks  to  elevator. 
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The  final  section.  Section  V,  surmarizes  the  major  identities  and  relation¬ 
ships  developed  in  this  report.  The  results  have  been  put  in  a  series  of  tables 
to  provide  a  ready  reference. 


SECTION  II 

MODAL  EESP0N3E  COEFFICIENTS  AND  SENSITIVITIES 

A.  MODAL  RESPONSE  COEFFICIENTS  AND  SYSTEM  RESPONSES 

For  the  single-loop  system  shown  in  Fig.  1  the  feedback  system  analysis 
process  of  interest  here  can  be  presumed  to  begin  with  a  given  nominal  open- 
loop  transfer  function  G(s)  in  factored  form.  The  results  of  the  next  stage  of 
analysis  are  the  closed- loop  transfer  functions  >  which  can 

also  be  in  factored  form.  The  closed-loop  system  has  but  one  input,  r(t),  and 
two  observable  outputs,  c(t)  and  e(t) ,  so  these  two  closed- loop  transfer  fvuac- 
tions  together  constitute  a  conplete  mathematical  model  for  the  system.  In 
essence  the  two  open- loop  elements  of  Fig.  2  replace  the  single  open-loop  trans¬ 
fer  function  G(s)  plus  the  feeaback  connection  (Fig.  1 )  as  a  system  model. 


System  Output 


System  Error 


Figure  2.  Open-Loop  Elements  for  Single-Loop  Feedback  System 

The  closed-loop  transfer  functions  can  be  combined  with  transforms  of  input 
functions  to  compute  output  time  responses  using  routine  inversion  procedures. 
The  system’s  natural  modes  are  explicitly  displayed  by  the  qj^'s  of  the  transfer 
function  denominator  factors,  but  the  relative  magnitudes  of  these  natural  modes 
are  indicated  only  implicitly.  A  basic  set  of  these  magnitudes,  or  modal 


response  coefficients^  are  explicit  featiires  of  the  weighting  functions,  which 
represent  the  system  time  responses  to  a  unit  iii5)ulse  input.  Thus,  a  complete 
physical  picture  of  system  response  characteristics  is  most  easily  achieved  when 
both  the  transfer  functions  and  weighting  functions  (or  their  constituent  modal 
.response  coefficients)  are  known.  To  form  a  complete  analysis  structure,  it 
should  he  possible  to  find  the  coefficients  from  either  the  open- loop  or  the 
closed-loop  transfer  functions,  i.e.,  they  should  be  exp].icitly  defined  by  the 
elements  contained  in  either  Fig.  1  or  2. 

In  the  present  article  modal  response  coefficients  will  be  related  to  open- 
and  closed- loop  transfer  functions  for  systems  containing  only  first-order 
closed-loop  poles.  Several  useful  connections  between  the  m+n  modal  response 
coefficients  will  also  be  developed.  The  coefficients  will  then  be  used  as 
elements  in  system  weighting  functions  to  show  their  connection  with  system 
responses  to  various  siii5)le  time  function  inputs.  Subsequent  articles  in  this 
section  will  relate  the  coefficients  to  sensitivity  measures,  and  will  generalize 
the  developments  to  systems  containing  multiple-order  poles.  The  actual  determi¬ 
nation  of  the  coefficients  from  the  several  forms  of  transfer  function  represen¬ 
tation  (e.g..  Bode  and  pole-zero  plots)  will  be  deferred  to  the  next  section. 


1 .  Modal  Response  Coefficients 

The  system  input-output  weighting  function  is  the  output,  c(t),  response 
when  the  input,  r(t),  is  a  unit  impulse.  Since  the  Laplace  transform  of  a  unit 
impulse  is  unity,  this  weighting  function,  w(t) ,  is  simply  the  inverse  Laplace 
transform  of  the  closed-loop  system  transfer  function  G2,(,(s), 


w(t)  =/''^[Grc(s^ 


(1) 


If  the  nvunber  of  poles  of  greater  than  the  number  of  zeros  (m  >  1 ) 
and  the  closed- loop  poles,  -q^^,  are  all  first  order,  then  a  partial  fraction 
expansion  of  Gy^(s)  results  in 


w(t) 


m+n 

L 


Li=l 


% 

s  +  q^ 


m+n 

=  E 

i=1 


(2) 
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is  the  ith  modal  response  coefficient.  It  may  he  evaluated  from 


Qi  =  j(s  +  q,i)Gj.(,(s^  ^ 


(3) 


If  m  =  0,  the  bracketed  quantity  of  Eq  2  has  the  constant  term  fc/l  +  k)  added, 
and  the  corresponding  time  response  contains  an  additional  term  K5(t)y^(l  +  k) . 
For  real  systems  m  >  1  always,  if  higher  frequency  characteristics  are  taken 
into  account;  but  low  frequency  approximations,  for  which  m  =  0,  are  often  used. 

As  an  aside  from  the  main  argument,  it  is  worth  noting  that 
splitting  Grc(s)  into  partial  fractions,  as  in  Eq  2,  is  equimlent  to 
replacing  the  (m  +  n)th  order  differential  equation  for  w(t)  by  m  +  n 
first-order  differential  equations  of  the  form 


Wi(t)  +  qiW^Ct)  =  Qi5(t),  i  =  1,  2,  ...,  m  +  n 


where  Wj^(O)  =  0  and  the  total  weighting  function  is 


w(t)  =  w^(t)  +  W2(t)  +  ..♦  +  Wi(t)  +  •••  Wj^+n(t)  (5) 


Thus  the  total  weighting  function  is  equivalent  to  the  summed  responses 
of  m  +  n  first-order  systems  excited  by  an  iiopulse  input,  as  shown  in 
Fig.  5*  The  modal  response  coefficients  can  also  be  thought  of  as 
initial  conditions  on  unity-gain  elemental  systems  (the  Qj^'s  in  Fig.  5 
replaced  by  1  ’s)  which,  with  no  other  excitation  (no  impulse  input  to 
the  system),  results  in  a  system  output  equal  to  the  weighting  function. 
For  this  interpretation  the  elemental  differential  equations  would  be 


Wi(t)  +  qiWi(t)  =  0,  i  =  1,  2,  •••,  m  +  n 


(6) 


where 


Wi(0)  =  % 


The  modal  response  coefficients  can  also  be  evaluated  directly  from  the  open 
loop  transfer  function.  In  terms  of  G(s),  Eq  5  can  be  written  as 


Qi  = 


Rs  +  q^)G(s) 
L  1  +  G 


G(s)  J 
\  -^s=-q^ 


(7) 
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Figure  3*  Elemental  First-Order  Systems  Corresponding  to 
Components  of  the  Weighting  Function 


At  this  point  it  is  convenient  to  introduce  a  new  variahle,  Q^(s),  which  is 
defined  as 

,  (s  +  qi)G(s) 

=  ■  i"+G(sl  ■ 


and  from  Eq  ^ 


=  % 


Rewriting  Eq  8  as  jj  +  G(s^  ^(s)  =  (s  +  qj^)G(s) 


and  differentiating  with  respect  to  s,  gives 


D  -  «(  =  )]  ^  <«(S)  ^  =  0(3)  .  (S  .  ^ 


Evaluating  Eq  10  at  s  =  -qj^  and  recalling  that  G(-q^)  =  -1 , 


Qi  =  ^(-qi)  = 


-1 

dG(s)/ds 


Besides  providing  a  means  for  evaluating  the  Q^^'s  directly  from  the  open- loop 
transfer  function,  Eq  11  will  subsequently  be  used  to  show  that  a  simple  rela¬ 
tionship  exists  between  the  modal  response  coefficients  and  the  gain  sensitivities 
defined  in  Article  B  of  this  section. 

The  modal  response  coefficients  for  the  input-output  weighting  function, 
w(t),  are  closely  connected  with  similar  coefficients  for  the  input-error  weight¬ 
ing  function,  Wj,g(t).  In  general  the  error  is  given  by 

e(t)  =  r(t'1  -  c(t)  (12) 


so  the  error  re^onse  to  a  unit  impulse  input  will  he 


Lvp 


Vj,e(t)  =  8(t)  -  w(t) 

1^+“  -a  t 
S(t)  -  E  ,  m  ^  1 

i=l 

^  -qit 

5(t)  -  ]C  QiS  ,  m  =  0 


(15) 


The  modal  response  coefficients  for  the  input-output  weighting  function  are  thus 
seen  to  be  the  negatives  of  similar  coefficients  for  the  input-error  weighting 
function. 


m 


»•••*  "jV 


J  -  i 

PS? 


Certain  simple  combinations  of  the  modal  response  coefficients  have  properties 
which  are  occasionally  useful.  These  are  most  conveniently  developed  by  matching 
coefficients  in  expansions  based  on  the  partial  fraction  form  and  the  closed- loop 
system  e:q)ressions  given  as  part  of  Fig.  1 .  In  root-locus  form,  when  m  >  1 , 


PC  n  (s  +  Zj) 

P  _  Kctjs)  _  j=1 _ 

rc'  '  ~  p(s)  +  Ka(s)  m+n 

n  (s  +  qi) 
i=1 


(1^^) 


The  sum  of  the  negatives  of  the  closed-loop  poles  is  (Ref.  1 ) ,  for  m  >  1 , 


lUTll 

£  q  =  "bi  + 
j=1 

where  is  the  sum  of  the  negatives  of  the  open- loop  poles 


(m+n  . 

”1  =  g 


and  6. 


is  the  Kronecker  delta  (6^  =  0  for  i  ^  i  =  j)*  Equation  l6  then 


becomes 


m+n  1  m+n  or  la+n  -I 

s"™-'  2  %  +  s"“-2r{b,  +  o  2  9i  -  E  ftiiil  +  •••  =  +  • 

i=l  •-  i=1  i=l  J 


Ey  matching  coefficient s. 


m+n 

E  %  =  '«m 

i=1 


m+n  m+n 

(b-j  +  k8iii)  ••  ^i^i  ~  +  iO.-\ 

i=1  i=1 


Inserting  Eq  I9  into  Eq  20 


m+n 

L  Qiii 
i=l 


=  K(b^  -  a^  +  /c)6ji  -  k6^ 


The  results  of  Eq  19  and  21  are  summarized  below. 


The  Bode  form  for  Gj^(s)  can  also  he  used  to  establish  similar  results. 
OSie  simplest  of  these  occurs  when  s  is  allowed  to  approach  zero,  i.e.,  for 
n  >  1, 


lim  Qj  l-im  KA.(s) 

=  4i  "  s>^B(s)  +  KA(s) 


(22) 


or 


m+n  Q. 

E  —  = 

i=i  % 


1 

K 


,  k  >  1 

1  +  K  ^  ^  ° 

0  ,  k  ^  -1 


(25). 


Equations  1 9,  21 ,  and  25  can  be  used  to  calculate  one  or  two  of  the  ' s  if  all 
the  others  are  known,  or  as  checks  on  computation. 

2.  System  Responses  to  Simple  Inputs 

Besides  their  central  role  in  the  system's  weighting  functions,  the  modal 
response  coefficients  also  appear  in  simple  ways  in  the  system  response  to 
other  elementary  inputs.  To  illustrate  this  feature  the  outputs  resulting  from 
steps,  periodic  inputs,  and  inputs  composed  of  power  series  in  time  will  now  be 
examined. 

S’*  Response  to  Step  Inputs.  Next  to  the  weighting  function  the  step 
response  is  the  most  common  transient  response  model  used  in  systems  analysis. 
The  indicial  response,  or  response  to  a  unit  step  input,  is  just  the  integral 
of  the  weighting  function. 


where  the  Qi^'s  are  the  partial  fraction  coefficients  for  the  step  response,  m> 
and  the  q- 's  are  all  first  order. 


The  Qj^'s  are  related  to  the  modal  response  coefficients,  Qi’s,  by 


^is  [s  1i)'^rc(^^ 

•-  -*  s=-q.i 


(25) 


The  amplitude  of  the  step  con^onent  in  the  system  oui5)ut,  Q^,  is 


1  ,  k  ^  1 


K 


%  -  [Grc(®)]  "  1  1  +  K  ^  ^  ° 


0  ,  k  <  -1 


or,  comparing  Eq.  26  with  Eq  23, 


(26) 


m+n  Q. 

=  Csr=(=)]  ^  ^ 

*-  s=0  i=1  % 


(27) 


Qq  is  the  magnitude  of  the  output’s  final  value  when  the  input  is  a  unit  step 
and  the  system  is  stable  (real  part  of  q^  >  O) .  Combining  Eq  24,  25,  and  27, 

m+n  Q.  -o.t 

I(t)  =  I)  ^  (1  -  e  0  (28) 

i=1  ^i 

b.  Response  to  Periodic  Inputs.  Any  periodic  input  may  be  represented  by 
Fourier  series  of  sine  and  cosine  terms  at  various  frequencies.  Because  the 
systems  considered  here  are  linear,  the  output  is  the  sum  of  the  outputs  due  to 
each  term  in  the  Fourier  series.  Consequently,  only  one  sine  and  one  cosine 
input  need  to  be  considered.  To  modify  the  periodic  input  (which  is  defined 
from  Si  t  ^  from  a  Fourier  series  standpoint)  to  a  periodic -like  input 
which  is  zero  until  t  =  0,  the  elemental  periodic  input  will  be  multiplied  by  a 
unit  step,  u(t) .  The  input  is  then 


(29) 


r(t)  =  cos  o^t  +  Tjj^  sin  a:^-^u(t) 
=  [cjj.  cos  (o^t  -  lifjj.)]  u(t) 


where  Cjj.  =  ^af  +  h|,  ilfjj.  =  tan“^ 

The  respons.e  for  an  arbitrary  input  is  related  to  the  impulse  response  by 
the  convolution  integral,  i.e., 

c(t)  =  I  r(t  -  T)w(T)dT 


-P 

-P 


r(T)w(t  -  T)dT 


With  the  weighting  function  expressed  in  terms  of  the  modal  response  coefficients 
(Eq  2)  and  the  input  given  by  Eq  29,  the  output  becomes,  using  Eq  30, 


c(t)  =  2^  - r - r -  e 

i=1  (if  + 


m+n 

£  [aj,  cos  (o^t  -  q^i)  +  bj^  sin  (o^t  -  (51) 

i=1  */qf  +  (xt 


1=1  yii  + 

where  <Pj^^  =  tan"^  ^/‘li*  terms  of  Cj^  and  the  steady-state  portion  of 
Eq  31  (the  transient  first  term  is  unchanged)  is 


m+n 

css(t)  =  ^k  D  -/- --  -  ^k  -  "Pki) 

i=1  yq|  +  ci^ 


The  ith  terms  in  the  responses  of  Eq  51  an<i  52  are  just  the  response  of  the  ith 
elemental  first-order  system  shown  in  Fig.  5  to  the  periodic-like  input  given 
by  Eq  29.  Considering  only  the  steady-state  portion,  each  elemental  first-order 
system  scales  the  input  amplitude  by  the  amplitude  ratio  Qj^/ ^q2  +  o^,  and  con¬ 
tributes  a  phase  lag  9^^^. 


.***  /- 


Equation  51  or  52  can  also  "be  used  to  prove  the  well-known  fact  that  the 
amplitude  ratio  and  phase  of  the  steady-state  response  to  a  sinusoidal  input  of 
frequency  is  the  anrolitude  and  phase  of  Gj.g(ja3i^)-  Thus,  in  more  conventional 
form. 


Css(i^)  =  Ci.lGrc(ja^)|  cos  [o^t  -  ilTj,  +4rGrc(J«k)]  (55) 


where 


m+n  0^ 

!G'vc(j^)l  cos  ^Gj.qv jo^)  — 


i=l  l^qf  +  4 


cosq)^^.  =  D 


and 


m+n 

lGvc(j"^)l  sin^G  (j%)  =  -  D  /  2  2  " 

i=i  yqf  +  4 


i=l  qf  + 


^  Qi% 


i=l  qf  +  4 


For  a  periodic-like  input  containing  all  the  terms  of  a  Fourier  series,  i.e. 


r(t)  =  [a^  cos  o^t  +  hj^.  sin  a:^^u(t)  (5^) 

k=0 


the  response  will  he,  by  superposition. 


^  -  ^k^i) 

c(t)  =  2^  D  { - r - 5 -  e 

k=0  i=1  ( 


% 


2 

+  0^ 


[fk  cos  («k^  -  'Pki)  \  I 


(35) 


c. 


Response  to  Power  Series  Inputs.  For  inputs  which  can  be  represented 


i{t)  -  %  S  M+1  r  j-J 

1=1  (-q,-)  *-  0=0*^-* 


Thus,  for  an  input  of 


r(t)  =  E  ^  ,  t  ^  0 

=  0  ,  t  <  0 


the  resDonse  is 


=  2:  2^  ^  ‘  ~1^~  ^ 

i=l  M  (-qj^)*  L  j=0  J*  J 


q^  r  0 


For  exanrole,  the  response  to  a  unit  step  (M  =  0)  is 


m+n  Q.  /  -q.t  \ 

=  ■’)  ’ 


which  corresponds  to  the  previous  results  (Eq  28) .  For  a  unit  ranro  input. 


[c(t)] 

•-  '  Unit  ramp 


li^n  0,  I  -q.t 


= 

i=i  qf  ' 


I  +  q. 


it'j  ,  q .  0 


The  modal  response  coefficients  can  he  related  to  conventional  output  and 

error  coefficients  hy  using  only  the  portion  of  Eq  38  v^hich  contains  the  povrer 

-qi  t 

series  terms  in  the  output,  i.e.,  hy  ignoring  the  e  transient  terms.  Then, 
for  a  power  series  input 

r(t)  =  dQ  +  d^t  +  ^  t^  +  ^  +  •••,  t>0  (4i; 


the  output,  less  transient  tenns,  will  he 


m 


m 


u- 


>v. 

-y 


LjL*>  • 

v'-!. 


m 


ll'-l 


iA 

r: 


[-.* 


[c(t)3 


Power  series  only 


A  ,  -  iMH  ^  i! 


i=1  M=0  (-q.i)  j=0  j- 


(42) 


By  changing  the  summation  index  j  to  k  =  M  -  j  and  rearrayiging  the  order  of  the 
summations^  Eq.  42  becomes 


[c(t3 


Power  series  only 


oo  m-Hi 

-E  E 


% 


E  H 


M-k 


k=0  i=1  M=k  (M  -  k)l 


(45) 


From  Eq  4l  it  can  be  seen  that 


;-k 


d^r(t) 


M=k  (M  -  k)l 


dr 


(44) 


Combining  Eq  43  and  44, 


[c(t)] 


Power  series  only 


'h,  W 


m+n 

=  r(t)  ^ 
i 


at-n  Q.  m+n  u.  mt-n  o- 

E  ^  -  Ht)  E  \  *  Ht)  E  i  + 

1=1  'k  i=i  li  i=i  qf 

=  EQr(t)  +  K^r(t)  +  K2f(t)  +  ••• 


m+n  Q. 


(45) 


The  R's  in  Eq  45  are  output  response  coefficients  which  are  occasionally  useful 
for  special  purposes.  Far  more  common  are  the  error  coefficients.  The  power 
series  portion  of  the  error  response  to  a  power  series  input  is 


[5(^3  Power  series  only  power  series  only 

=  EQr(t)  +  E^  r(t)  +  E2r(t)  +  • • •  (46) 


21 


'-V- 

-V- 

Cv- 


I 

'.<K- 


j"' 


t 


fiV, 

Sesjs 


i'  W* 


I-*.'-* 


y.->. 


H-'V' 


K/: 


yvt'' 


;v 


•  *.  ** 

-  .  >  *-> 


Using  the  results  of  Eq  45,  the  error  coefficients  are  seen  to  he 


m+n  n 

Eq  =  1  -  Bo  =  1  -  E  - 

1=1  ^1 


E,  =  -R,  = 


m+n  Q. 

E  ^ 

i=l  *3.1 


Eg  =  -Rg  = 


m+n  0^ 

~5 
i=i  q? 


(47) 


%  =  -%  =  (-1) 


,  ,,  m+n  n. 
k+1  ^  ^ 

1  k+1 

1=1  q^ 


Because  of  the  many  ways  in  which  modal  response  coefficients  can  be  obtained 
these  formulas  are  especially  handy  for  the  calculation  of  error  coefficients. 

B.  SENSITIVITY 

The  modal  response  coefficients  discussed  above  are  the  last  step  in  the 
analysis  problem  for  linear  feedback  systems  if  the  open-loop  transfer  function, 
G(s),  is  considered  to  be  an  exactly  specified  quantity.  Unfortunately,  the 
poles,  zeros,  and  gain  of  G(s)  ordinarily  vary  somewhat  about  nominal  values 
which  are  themselves  inexact.  Thus,  analyses  are  not  really  complete  without 
some  consideration  of  the  effects  of  parameter  variations.  This  purpose  is 
served  by  the  sensitivity  factors  to  be  discussed  in  this  article. 

As  noted  in  the  Introduction,  sensitivity  measures  have  been  a  major  concern 
in  feedback  systems  since  the  early  days  of  feedback  amplifiers.  The  pioneering 
investigators  defined  a  sensitivity  function  as  one  of  the  essential  features  in 
a  mathematical  definition  of  feedback.  This  definition,  which  relates  to  the 
over-all  sensitivity  of  closed-loop  transfer  functions  to  variations  in  the  open- 
loop  transfer  function,  we  shall  call  "classical  sensitivity."  More  recently 
attention  has  been  centered  on  the  effects  upon  closed- loop  poles  of  variations 
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in  open-loop  gain,  pole,  and  zero  locations.  The  sensitivities  of  this  nature 
we  shall  call  "gain  sensitivity,”  "pole  sensitivity,"  and  "zero  sensitivity," 
respectively.  These  sensitivities  are  especially  important  in  servosystems 
because  most  such  systems  exhibit  only  one  or  two  dominant  closed- loop  modes. 

The  variations  of  these  modes  with  open-loop  characteristics  are  therefore  of 
most  concern,  and  the  pole,  zero,  and  gain  sensitivities  are  of  most  value  in 
assessing  the  effects  of  the  variations.  So,  the  sensitivity  functions  for 
closed-loop  root  variations  will  be  the  major  topic  below,  with  their  relation 
to  classical  sensitivity  beirig  a  short  aside.  Again  the  developments  will 
presume  first-order  closed-loop  poles,  with  generalization  to  the  multiple-order 
sitviation  being  deferred  to  the  last  article  of  the  section. 


1 .  Gain,  Pole,  and  Zero  Sensitivities 


The  purpose  of  this  article  is  to  derive  the  expressions  for  the  first-order 
sensitivity  of  the  closed- loop  poles  to  variations  in  the  open-loop  gain,  fc,  the 
open- loop  poles,  -p.,  and  the  open- loop  zeros,  -  ..  VJith  all  these  parameters 
being  considered  as  variable,  the  open-loop  transfer  function,  G,  must  be  con¬ 
sidered  as  a  function  of  all  of  them,  i.e., 

G  =  G(s,fC,Zj,Pj)  (48) 

Then  the  total  differential  of  G  is 


8g  ,  8G  ,  ^  3G  ,  8G  , 

,1=1  j  j— 1  J 


(49) 


Because  the  closed- loop  poles,  -q^,  are  defined  by  the  equation 


D  *  °(^)]  ° 


(50) 


the  total  differential  of  Eq  48  must  be  zero  for  s  =  -q^.  Setting  dG  =  0  and 
s  =  -qj^  in  Eq  49,  and  rearranging  terms  gives 


S^‘ir.4jd 


=  [5^]  |(i)  (^)  ^  <5l) 

*-  '  -'s=-qj^l  s=-qj|^  j=1  '  0's=-q^  j=l  '  ^J's=-qjj^  ) 


!rhe  variation  in  the  negative  of  the  closed- loop  pole,  dqj^,  can  also  he 
expressed  in  anocher  way  hy  noting  that  q^  depends  only  on  K,  the  Zj*s,  and  the 


Pj  i*e.^ 


li  = 


and  writing  another  total  differential  as 


dq.  = 


^^i  dK  ^  ^  °^i  ,  ^ 

Ss: 


^  ^1i 


m+n  Sq. 


dp. 


.  n  .  m+n  . 

=  —  +  ^  S2  dz-  +  si  dp 

K  K  ~  2.  j  p. 


<5=1 


j=1 


(52) 


Here  the  factors  Sq^^/Bz^,  and  ^q^^/^Pj  are  given  a  special  symbol,  S. 

These  are  the  first-order  sensitivity  factors.  The  subscript  and  superscript 
notation  indicates  that  a  differential  increment  in  the  open- loop  quantity, 
defined  by  the  subscript,  resijilts  in  a  differential  increment  of  the  ith  closed- 
loop  root  equal  to  the  sensitivity  factor  times  the  open-loop  parametric 
variation. 

Equating  like  coefficients  in  Eq  51  and  52  gives 


SS-i«R» 


K 


Sqi 

^k/k 


=  K 


^G/ds 


s— q. 


foG/dz^ 

\  dG/6s  I 


dq^  / dG/Sz . 

Sqi 


(53) 


( ^G/apj^ 


\  dG/bs 


s=-q. 


s=-q. 


Note  that  the  gain  sensitivity  is  based  on  a  fractional  change  in  K,  while  the 
pole  and  zero  sensitivities  are  based  on  absolute  shifts  of  pj  and  zy  These 
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definitions  were  selected  because  they  provide  some  simplifications  in  subse¬ 
quent  relationships. 


If  the  open- loop  transfer  function,  G,  is  written  in  root-locus  form. 


n  (s  +  Zj) 

G  =  K  - 

m+n 


n  (s  +  p.) 

j=l 


the  sensitivities  are  (remembering  that  at  s  =  -q^^,  G  =  -1 ) 
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Examination  and  further  interpretation  of  Eq  55  reveal  four  very  interesting 
properties  of  the  sensitivity  factors; 

a.  The  gain  sensitivity  is  a  factor  in  each  of  the  sensitivity 
terms.  Thus,  Eq  52  becomes 

•  C  „  Q  dz .  m+n  dp .  ”1 

L  0=1  0  0=1  ^  oj 

b.  The  gain  sensitivity  is  equal  to  the  modal  response  coeffi¬ 
cient.  This  is  easily  seen  by  comparing  the  expression  for 
the  gain  sensitivity  with  Eq  11,  i.e.. 


%  =  [dG(s)/ds]^^_^^  = 


This  equality  is  very  useful  since  all  of  the  properties 
previously  derived  for  the  modal  response  cofficients  are 
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app3.icat)le  to  the  gain  sensitivity.  Using  this  correspond¬ 
ence,  other  formulas  for  the  gain  sensitivity  include 


Rs  +  ( 

"  L  1  + 


q.j,)G{s)1 


/C  n  (-qi  +  Zj) 

_  j=^ _ 

m+n 

(i  +  K&^)  n  (-qi  +  qj) 
j=i 

m+n 

-  n  (-qi  +  Pj) 

_ j£! _ 

m+n 

(1  +  Kbjn)  n  (-qi  +  qj) 

j=1 

Another  formula  for  S^,  which  could  have  been  derived  previ 
ously  for  Qji,  is  a  direct  result  of  Eq  57: 


Various  gain  sensitivity  combinations  have  the  same  simple 
forms  as  the  modal  response  coefficients.  These  follow  from 
the  equality  of  Eq  57  and  the  previous  results  of  Eq  1 9^  21 , 
and  25,  i.e.,  for  m>  1 


a 

i 


i 

I 

I 

/P 

'01 

i 


1 


I 


■~'l 


’• 

bi 

il:‘ 


rv 

I 

I 

XI 


d.  The  sum  of  all  the  zero  and  pole  sensitivities  for  each 
closed-loop  pole  must  equal  one.  This  follows  directly 
from  Eq  55  and  59*  Thus, 


n 

E  sL 

j=i  J 


m+n 

E  s 

ci=i 
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This  initially  surprising  result  is  easily  explained  by 
recalling  a  root-locus  plot.  If  all  the  open- loop  zeros 
and  poles  are  moved  the  same  amount,  all  the  closed- loop 
poles  will  be  moved  by  that  amount. 


(6i) 


The  gain  sensitivity  can  be  interpreted  in  two  ways  which  are  physically 
enlightening.  First,  is  a  measure  of  the  slope  of  a  conventional  root  locus. 
As  root  loci  are  ordinarily  plotted  for  fixed  open-loop  poles  and  zeros,  the 
only  variable  along  the  loci  is  gain.  Thus  the  plot  gives  the  zeros  of 


1  +  G(s,k)  =  0 


Taking  the  total  derivative 


dG 


^G 


ds  + 


hG 


dfc 


Along  the  locus  dG(s)  =  0,  so 


ds 


J 


s= 


d/C 

K 


(62) 

(65) 


(64) 


d/c//c  is  a  real  number,  so  the  direction  of  ds  along  the  locus  for  positive  d/c//c 
will  be  given  by  -S^  which,  in  general,  will  be  complex.  The  minus  sign  appears 
because  is  the  sensitivity  of  q^  which  is  the  negative  of  the  ith  closed- loop 
pole. 

The  second  interpretation  of  derives  from  the  same  argument  as  that  used 
above  by  adding  an  additional  step.  The  derivative  of  In  G(s)  is 


d  In  G(s)  ^  1  dG(s) 

ds  G(s)  ds 
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wMch,  along  the  root  locus,  hecomes 

d  In  G(s)  ^  _1_ 
ds  qi 

or 

oi  _  ds 
K  -  d  In  G(s) 

- 

~  dC/G 


(66) 


(67) 


Thus  the  gain  sensitivity  is  a  measure  of  the  shift  in  a  closed- loop  pole  due 
to  a  fractional  change  in  the  open-loop  transfer  function.  Because  only  K 
changes  along  a  conventional  root  locus,  Eq.  67  is  really  no  different  from 
Eq  6i^,  except  for  the  introduction  of  the  logarithmic  form. 

2.  Limiting  Behavior  and  Special  Cases 


The  magnitudes  of  gain  sensitivities  can  cover  the  entire  range  of  values 
from  minus  to  plus  infinity.  Yet,  intuitive  notions  of  "sensitivity"  as  a 
general  concept  in  closed- loop  systems  make  part  of  this  range  appear  unreason¬ 
able.  One  part  of  the  problem  is  a  direct  consequence  of  the  sensitivity 
definition,  while  another  is  associated  with  its  first-order  approximation  nature. 
A  better  vinderstanding  of  both  facets  can  be  gained  by  an  examination  of  limiti ng 
cases. 

In  general,  closed-loop  poles  depart  from  open-loop  poles  for  low  values  of 
gain,  and  proceed  to  either  open-loop  zeros  or  unbounded  values  as  the  open-loop 
gain  becomes  very  large.  The  gain  sensitivity,  as  given  by  Eq  59^  is 
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0=1 
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m+n 
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j=i 


(59) 


As  K  approaches  zero,  the  closed-loop  root  -q^  approaches  the  open- loop  pole 
from  which  it  derives,  i.e.,  q^^  ^  Pj^*  Then  the  term  l/(qjL  -  p^)  in  Eq  59  is 
dominant,  so 
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Similarly,  as  K  "becomes  very  large,  n  of  the  closed-loop  poles  approach  open- 
loop  zeros.  If  the  ith  closed-loop  pole  is  one  of  these,  and  it  approaches  the 
jth  open-loop  zero  so  that  then 


Jk-^c 


(^J  -  ^i) 


Finally,  m  of  the  closed- loop  poles  have  no  zeros  to  go  to,  and  hence  become 
very  large  relative  to  the  p.  and  z..  !Ehe  sensitivity  for  these  poles  is 

J  u 


ki  »  z,,p. 


“  1 


When  the  gain  is  sufficiently  large  for  the  open- loop  zero  db  line  to  intersect 
the  high  frequency  asyntptote,  the  open- loop  transfer  function  is  approximately 


G(s)  4  4 


so  that  qj^  will  be 


Thus,  the  sensitivity  of  the  unbounded  pole  will  be 


mr— 


Equation  72  indicates  that  the  sensitivity  increases  as  the  mth  root  of  K  as 
gain  is  increased,  although  for  finite  gains  the  sensitivity  is  always  finite. 


Another  circuEStance  in  vbich  the  sensitivity  can  oecosse  very  large  is 
revealed  oy  Bq  58d.  Here  it  is  apparent  that  the  gain  sensitivity  for  a  closed- 
loop  pole  becoces  very  large  as  that  pole  nears  another  closed- loop  pole. 

Indeed,  as  beccres  equal  to  q^,  indicating  a  branch  point  on  the  root  locus, 
the  gain  sensitivity  goes  to  infinity.  This  is  to  be  expected  since  the  sensi¬ 
tivity  factors  defined  thus  far  have  not  considered  multiple-order,  closed- loop 
roots.  As  long  as  the  gain  is  finite,  an  infinite  gain  sensitivity  always 
indicates  nultiple-order,  closed- loop  poles - 

A  special  situation  of  considerable  interest  can  occur  when  a  closed- loop 
root  lies  between  an  open-loop  pole  and  zero  which  are  much  closer  to  each  other 
than  to  all  other  open-loop  poles  and  zeros.  This  is  the  so-called  dipole  case. 
Sje  sensitivity  for  the  bounded  closed-loop  pole  will  be,  approximately. 


5-  Sensitivity  Functions  for  Alternate  Transfer  Function  Forms 

The  above  equations  for  gain  sensitivity  are  the  same  whether  G  is  written 
in  root-locus  or  Bode  form,  i.e.,  =  S^.  In  fact,  the  equations  are  still 

valid  if  some  of  the  terms  in  G  are  in  root-locus  form  and  some  are  in  Bode 
form.  It  is,  however,  necessary  to  modify  the  open-loop  pole  and  zero  sensi¬ 
tivities  for  terms  which  are  written  in  Bode  form.  From  Eq  55  it  can  be  shown 
that  for  zeros  and  poles  which  appear  in  G  in  Bode  form,  i.e.,  [(s/zj)  +  l]  or 
ps/pj)  +  the  sensitivities  are 


(75) 


^3  PjUi  -  Pj; 


Frequently  open- loop  zeros  or  poles  will  occur  as  con^lex-con jugate  pairs, 
and  variations  in  the  system  will  change  both  zeros  or  poles.  Consequently  it 
becomes  desirable  to  introduce  sensitivities  for  the  parameters  which  define  a 
complex  pair  of  zeros  or  poles.  For  example,  consider  a  complex  pair  of  zeros, 
z-j  and  Z2,  which  are  defined  by  their  frequency,  oi,  and  damping  ratio,  i.e.. 


z-|  =  ^0)  +  jcD  yi  - 

Z2  =  V’’  ~ 


For  this  situation  it  is  useful  to  define  frequency  and  dauqping  ratio  sensi¬ 
tivities  as 


,i 

Sqi  dz^  ^q^  dZ2 

bzj  do  'bz2  do 

,i 

Sq^  dz^  Sq^  dZg 

II 

dz-j  d^  ^  dz2  d^ 

It  is  easily  shown  from  Eq  55^  T6,  and  77  that  if  the  term  appears  in  G  in 
root-locus  form  (s^  +  2^a)s  +  oc?) ,  the  frequency  and  damping-ratio  sensitivities 
for  a  complex  pair  of  zeros  or  poles  are 


±2((o  -  ^q.  )sj 

qS  -  21^0X1^  +  ^ 

q.f  -  2^0X1^  +  0^ 


where  the  upper  sign  is  to  be  used  for  zeros  and  the  lower  one  for  poles. 


If  the  term  is  written  in  Bode  form  [(s2/a?)  +  (Sis/ai)  *  l], 
are 


the  sensitivities 


±2q^  (^co  - 


^  “  qf  -  ^ 


(79) 


Sj  = 


+2qj^ajSj^ 


^  qf  “  2^<uq.jL  + 


For  some  cases  it  may  he  more  convenient  to  define  a  complex  pair  in  terms 
of  their  real  and  imaginary  parts,  i.e.. 


2^  =  a  +  Jb 
22  =  a  -  jb 


In  this  case,  with  the  term  in  root-locus  form  +  2as  +  (a^  +  ,  the  sensi¬ 

tivities  are 
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If  the  term  is  in  Bode  form  |s^/(a^  +  b^)  +  2as/(a^  +  b^)  +  1 J ,  the  sensitivities 


are 
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b2) 
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)S^ 
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52 

(80) 


(81) 


4.  Connection  Between  Classical  and  Pole,  Zero,  and  Gain  Sensitivities 


As  pointed  out  by  Bode,  the  usual  conception  of  feedback  involves  two  distinct 
ideas.  The  first  is  the  more  obvious  and  connnon — that  feedback  implies  a  loop 
transmission  or  return  of  some  measure  of  output  quantities  to  earlier  stages  of 
the  system.  The  second  notion  is  that  of  a  reduction  in  the  effects  of  varia¬ 
tions  in  the  elements  of  the  forward  loop. 

To  point  up  these  old  and  well-known  singale  ideas,  consider  the  system  shown 
in  Fig.  4.  Here  Gg^  is  the  main  controller,  Gp  the  feedback  element,  Y  the 
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output,  C(s),  is 


C(s)  =  — 


^a^  Y 

■  GaGfY  ^  1  +  GaGpY 


=  G,c[E+i»] 


=  Gi,e|GaYR  +  YN] 
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1  +  GaGfY  1  +  G 

The  presence  of  the  loop  transmission,  i.e.,  the  existence  of  the  feedback 
element  G-p,  can  lead  to  practical  advantages  when  feedback  is  dominant.  Thus, 


wiien  |gJ  1  y  the  outout  approxiiaates 
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ae  relationship  heteeen  output  and  input  is  rade  substuntialiy 
of  the  eontroUed-eleinent,  and  the  effect  of  unvanted  signals  on  tee’output  is 
mtenally  reduced  (%  being  generally  very  large  to  nake  possible  the  inequality 
^ •  As  an  open-loop  control,  the  output  would  be 


C(s)  =  (^YR  -f  YIJ 


rne  difference  due  to  the  closing  of  the  loop  is  just  the  error/input  tr^sfer 
function  C-re-  Thus  Gj-e  (the  classical  "return  difference”)  is  the  fuitoental 
measure  of  the  iinprovement  introduced  by  the  loop  tiansmission. 


To  illustrate  the  decreased  sensitivity  effects,  consider  a  differential 
change  in  the  closed-loop  transfer  function  G^c.  For  regions  in  the  s-p3ane 
where  Gj.g  and  G  are  analytic. 
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or,  in  terms  of  fractional  changes. 


=  1  ^ 

“  1  G  [Ga  Y  '  ^ 


LGa  Y  Gf  J 


Again,  when  the  feedback  effect  is  dominant. 


^^rc  .  T_  ^  dGf 

Grc  G  [_Ga_  Yj  Gf 


The  beneficial  consequences  of  large  values  for  the  open-loop  transfer  function, 
G,  in  decreasing  the  effect  on  the  over-all  transfer  function  of  forifara-laop 
controller  and  controUed-elensnt  variations  is  ap^rent  flora  Bq  67-  Sis'  open- 
loop  variatioisdC^(^  and  dY/Y  are  decreased  by  the  factor  i/G.  Kote,  hor-ever, 
that  variations  in  the  feedback  eleraents  are  reflected  directly  into  variations 
of  the  closed-loop  transfer  function  Gj^-  In  general,  the  beneficial  reduction 
is  deterained  by  the  factor  l/(l  -r  G) — again  the  error  transfer  function  is  a 
fusdanental  ceasure,  this  tine  of  sensitivity. 

itefine,  now,  an  over-all  systen  sensit.i.vity  as 


J^xc  ^^rc/^rc 


d  In  G^ 
d  In  X 


^  ^  Grc 
dx/x 


(88) 


where  x  is  soise  caiarseter  or  element  in  the  open-loop  systen.  If  x  is  taken  to 
be  the  forward-loop  controller  transfer  function,  G^,  the  controUed-elenent 
transfer  function,  Y,  or  the  forward-loop  gain,  Xp,  it  is  apparent  that 


This  definition  of  sensitivity  may  be  considered  to  be  classical  (Ref.  p  a-'*-’! 
although  Bode’s  original  definition  amounted  to  the  inverse  of  that  defined  here. 
Examination  of  the  role  played  by  the  classical  sensitivity  in  Eq  82,  86,  and  89 
indicates  the  overwhelming  importance  of  classical  sensitivity  as  a  fundamental 
quantity  in  feedback  systems. 

Classical  sensitivity  can  also  be  used  in  expressions  involving  open-loop 
gain,  pole,  and  zero  changes.  To  illustrate  this,  return  to  the  simplified 
block  diagram  of  Fig.  1  for  which 
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'G  "  dG/G  "  1  +  G  “ 


(90) 


The  fractional  change  in  the  open- loop  transfer  function  is  given  by 
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so  the  fractional  change  in  the  closed-loop  system  characteristics  is  just 
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=  a 
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(92) 


!nie  classical  sensitivity  is  a  function  of  s,  and  can  range  over  extreme 
values-  In  the  frequency  wosain  (s  =  jto)  its  lira  ting  cases  occur  when  JgJ  »  1, 
wiiich  usually  occurs  for  lo'w  frequencies,  and  |c-|  «  1 ,  which  always  occurs  for 
high  frequencies  (n  >  l),  viz: 
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(G(j(D)j 
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|G(ja>)|  »  1 
|G(ja))|  «  1 


(93) 


Therefore,  on  a  classical  sensitivity  basis  a  "good"  system  is  one  with  lots  of 
feedback — |G(ja))|  »  1  over  a  broad  frequency  range — ^hardly  a  revelation  of 

earth-shattering  originality I 


The  classical  and  gain  sensitivities  are  related  in  a  very  simple  way.  The 
relationship  is  most  easily  developed  by  noting  that  classical  sensitivity 
corresponds  to  the  input-error  transfer  function  which,  in  turn,  can  be  written 
in  terms  of  the  modal  response  coefficients. 
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Since  the  gain  sensitivities  are  equal  to  the  modal  response  coefficients. 
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Equation  95  expresses  the  classical  sensitivity  as  a  weighted  sum  of  gain  sensi¬ 
tivities.  This  could  he  carried  further  formally  be  introducing  pole  and  zero 
sensitivities  to  replace  the  one  (see  Eq  6l),  but  this  seems  hardly  worthwhile. 
Equation  9p  also  indicates  that  the  gain  sensitivities  are  the  negatives  of  the 
residues,  evaluated  at  the  closed- loop  poles,  of  the  classical  sensitivity 
(Ref.  10  and  12). 


C.  EFFECTS  OF  Kth-ORDER  CLOSED-LOOP  POLES 

The  developments  of  Articles  A  and  3  are  restricted  to  cases  for  first-order 
closed- loop  poles.  This  article  considers  the  additional  con^lexity  introduced 
by  the  existence  of  Nth-order  closed-loop  poles.  VJith  the  exception  of  N  =  2, 
Nth-order  closea-loop  poles  seldom  occur  in  practice.  The  developments  of  this 
article  are,  therefore,  somewhat  academic,  although  required  for  completeness. 
This  article  is  basically  indepenaent  of  the  main  body  of  the  report,  so  the 
casual  reader  my  skip  it  without  seriously  detracting  from  his  understanding 
of  other  material.  The  procedure  in  this  article  will  be  to  retrace  most  of 
the  prior  developments,  noting  what  additions  or  modifications  are  required  to 
generalize  the  results  in  the  case  of  an  Nth-order  closed- loop  pole. 

If  -qj^  is  an  Nth-order  pole,  the  partial  fraction  expansion  of  the  impulse 
response  or  weighting  function,  Eq  2,  will  contain  the  terms 


where  the  N  modal  response  coefficients  are  evaluated  by 


There  is  a  useful  relationship  between  Qijj  and  a  modified  gain  sensitivity, 
which  will  be  developed  shortly.  To  show  the  relationship,  it  is  convenient  to 
introduce  the  variable 


where 


Rewriting  Eq  98  as 


^  ■  1  +  G(s) 


ei„(-ii)  =  % 


[1  +  G(s)]^  (s)  =  (s  +  <lj)’*G(s) 


and  differentiating  with  respect  to  s,  gives 


[1  +  G(s)] 


^  *  qi)"‘'0(s)  +  (s  +  q,)“  ^  (100) 


Evaluating  Eq  100  at  s  =  -qj^. 


therefore 
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Repeated  differentiation  of  Eq  1 00  shows  that 
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=  0  ,  1  <  k  <  N-1 
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(105) 
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s=-q. 
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Equation  105  expresses  the  interesting  property  that  when  -q^  is  an  Nth-order 
pole,  the  first  (N  -  1 )  derivatives  of  G(s)  with  respect  to  s  are  zero  at 
s  =  -q^.  Equation  104  is  the  general  expression  for  Qijj  and  reduces  to  Eq  11 
for  N  =  1 .  This  equation  will  be  used  to  show  the  relationship  of  Qijj  to  the 
modified  gain  sensitivity. 

A  relationship  between  the  N  modal  response  coefficients  can  be  derived 
from  the  above  by  considering  the  functions 


where 


(105) 


Using  an  alternative  expression  for  the  closed- loop  transfer  function  (Fig.  l), 


Eq  105  can  be  written 
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where 


m+n 
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(S  +  qj) 


product  of  s  minus  each  closed-loop  pole  except  the 
Nth-order  one,  -q^^  (there  are  m+n-N  terms  in  the 
product) 


After  differentiating  once,  Eq  106  becomes 
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Combining  Eq  107  with  the  formula  for  the  derivative  of  a  product^ 
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Combining  Eq  106  and  109^ 
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Once  Qi^  has  been  evaluated  by  one  of  the  methods  to  be  presented  in  Section  III 
Eq  111  can  be  used  to  calculate  the  remaining  (K  -  1 )  coefficients  for  the  Nth- 
order  pole,  -q^.  Note  that  the  coefficients  must  be  computed  in  the  sequence 


Qi  j  Qt  }  Qi  •  •  •  Qi ,  •  A  somewhat  similar  relationship  is  given  in 
^N-l'  ^N-2"  ^K-5 


Ref.  llf. 


Although  the  general  expression,  Eq  1 1 1 ,  is  rather  complicated,  for  the  most 
common  case  of  N  =  2  it  reduces  to  simply 


n  m+n. 
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(112) 


For  Nth-order  closed- loop  poles,  the  previous  expression  for  the  sum  of  the 
modal  response  coefficients,  Eq  19^  is  readily  modified  to 
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(115) 


where  the  summations  include  the  Qj^’s  for  all  first-order  poles  and  the  's 
for  all  higher  order  poles.  The  validity  of  Eq  113  can  he  shown  very  easily 
from  the  partial  fraction  expansion  of  Orc(s).  Equation  1 1 5  is  obtained  by 
multiplying  Gi.c(s)  and  its  partial  fraction  expansion  by  s  rnd  letting  s  -^oo. 


For.  Nth-order  poles  the  system  responses  are  considerably  more  complicated 
than  those  given  previously  for  N  =  1 .  For  example,  if  -q^^  is  a  second-order 
pole  the  system  response  to  an  input  of  sin  cut  contains  the  terms 
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For  an  input  of  cos  cut,  the  response  contains  the  terms 


r  1 

^2  ^  ^  jj^i®  ^  *^i  tut  +  o)  sin  oj^ 

" l4^ 

+  (q|  -  a^)  cos  ost  +  2qj^a5  sin  oitj 

A  general  expression  for  the  response  of  an  Nth-order  pole  to  a  cos  ojt  or 
sin  (ut  input  can  he  derived  by  considering  cos  cut  and  sin  cut  to  be  the  real  and 
imaginary  parts  of  e^“^,  iie.. 


cos  cut  =  Re  (e^^) 
sin  cut  =  Im  (e*^^) 


(11 4) 


Using  this  technique  we  find  that  response  of  an  Nth-order  pole,  -q^,  to 
cos  cut  input  is 
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For  an  input  of  sin  cut,  the  response  is 
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It  IS  important  to  note  that  in  evaluating  the  real  and  imaginary  parts  of  the 
above  expressions  q^  must  be  treated  as  a  purely  real  quantity  even  if  it  is 
actually  complex. 

The  responses  for  power  series  inputs  are  equally  as  complicated.  For 
example,  if  -q^  is  a  second-order  pole  other  than  zero,  the  response  to  a  unit 


step  will  contain  the  terms 


For  a  vinit  ramp,  the  terms  are 


M 

The  general  expression  for  the  response  of  an  Nth-order  pole  to  an  input  of  t  /M! 
can  be  derived  by  taking  a  partial  fraction  expansion  of  the  product  of  the 
Laplace  transform  of  the  input  and  the  terms  of  the  transfer  function.  After 
a  considerable  amount  of  manipulation,  the  inverse  transform  gives  a  response  of 
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From  the  previous  development  it  is  obvious  that  the  sensitivity  definitions 
must  be  modified  for  Nth-order  closed-loop  poles.  With  the  earlier  definitions, 
the  gain  sensitivity  is  equal  to  -(Bg/^s)"]^  ,  but  from  Eq  10^ 

s--qj^ 

°  ("5) 


or  the  original  gain  sensitivity  becomes  infinite.  To  modify  the  sensitivity 
definitions  for  Nth-order  poles,  we  expand  the  total  differential  of  G  to  include 
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Retaining  only  the  lowest  order  terms  for  each  parameter  and  remembering  that 
the  first  (ZJ  -  l)  derivatives  of  G  with  respect  to  s  are  zero  at  s  =  gives 
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s=-q. 


Equation  1 1 7  suggests  using  modified  sensitivity  definitions  such  that 
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From  Eq  117  and  II8  it  can  be  seen  that  if  G  is  in  root-locus,  Bode,  or  mixed 
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For  zeros  or  poles  which  are  written  in  root-locus  form,  the  sensitivities  are 
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For  zeros  or  poles  which  are  written  in  Bode  form,  the  sensitivities  are 
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Note  that  the  relationship  of  the  zero  and  pole  sensitivities  to  the  gain  sensi¬ 
tivity  is  exactly  the  same  as  that  for  first-order  closed-loop  poles;  Eq.  120  is 
identical  to  Eq  55 j  an*!  EQ.  121  is  identical  to  Eq  75*  I't  is  also  important  to 
note  that  from  Eq  1 04-  and  1 1 9 


(122) 


As  pointed  out  by  Ur  in  a  similar  development  (Ref.  15)»  Eq  11 8  reflects 
the  well-known  characteristics  of  a  branch  point  on  a  root-locus  plot.  The 
incoming  branches,  which  at  their  junction  represent  an  Nth-order  closed- loop 
pole,  are  evenly  spaced  and  separated  from  each  other  by  2jt/N.  The  outgoing 
branches  are  also  separated  from  each  other  by  2rt/N  and  are  midway  between  the 
incoming  branches. 


SECTION  III 


MODAL  RESPONSE  COEFPICIENr  AND 
GAIN  SENSITIVITY  EVALUATION  FOR  SINGLE-LOOP  SYSTEMS 

Modal  response  coefficients  and  gain  sensitivities  can  "be  obtained  in  a  wide 
variety  of  ways  from  either  open-  or  closed-loop  transfer  function  representa¬ 
tions.  The  representations  most  commonly  available  in  a  typical  problem  are: 

Open-loop  Closed-loop 

G(s)  in  factored  form 

Pole-zero  plot  Root  loci,  without  complete  sets 

of  compatible  roots 

G(ja))  Bode  plot  '»  Closed- loop  jUrcldb  l^rei^i) 

L  asymptotic  plots  (when  decompo- 

G(±ct)  Bode  plot  (also  serves  as  sition  method  is  applicable) 

a  plot  of  closed- 
loop  real  roots 
versus  gain) 

Correlating  the  infomation  appearing  on  all  these  common  representations  serves 
to  supplement  the  information  available  on  any  one.  For  example,  a  closed- loop 
|Gj.^,|db  asymptotic  plot,  found  using  the  G(±a)  and  G(jcD)  open- loop  Bode  plots 
and  the  decomposition  method,  yields  a  set  of  roots  which  can  be  noted  on  the 
root  locus.  Similarly,  one  branch  of  the  root  locus  complete  with  gains  might 
provide  sufficient  information  to  allow  a  decomposition  process  to  proceed  using 
the  Bode  representation.  In  any  event,  for  most  cases  in  practice  the  above 
forms  will  be  sufficient  to  provide  complete  information  about  gains  and  closed- 
loop  poles.  For  the  unusual  circumstance  where  this  is  not  true,  the  following 
additional  representations  may  be  necessary; 

Open- loop  Closed- loop 

G(s)  in  1  polynomial  form  Root  loci,  with  complete  sets 

of  compatible  roots 

G(5,n)  Bode  plots  Root  loci  as  functions  of  gain 

with  I  as  parameter 

These  various  forms  of  transfer  function  representation  supply  the  raw  data  from 
which  modal  response  coefficients  and  sensitivities  are  to  be  found. 


p*;#) 


The  following  discussion  of  modal  response  coefficient  and  gain  sensitivity 
evaluation  is  organized  along  traditional  lines,  with  an  individual  article  for 
each  general  form  of  transfer  function  representation.  The  articles  appear 
under  the  headings: 

A.  Direct  Calculation 

B.  Root-Locus  Methods 

C.  Methods  Using  Open-Loop  Bode  and  5  Plots 

D.  Method  Using  Closed-Loop  Bode  Asymptotes 

Ordinarily  different  formulas  for  modal  response  coefficients  or  gain  sensitiv¬ 
ities  are  most  appropriate  to  a  particular  form  of  representation.  However, 
because  of  the  supplementary  character  of  the  several  transfer  function  repre¬ 
sentations  available,  certain  formulas  work  well  with  several  representations. 
Consequently  some  interplay  at  the  detailed  formula  level  is  inevitable.  In 
practice  this  is  even  more  prevalent,  since  each  of  the  transfer  function  repre¬ 
sentations  available  in  a  given  problem  is  likely  to  be  most  suitable  for 
computation  of  a  particular  modal  response  coefficient.  Thus,  in  a  practical 
problem  one  might  get  several  coefficients  from  G(jto)  and  G(±a)  Bode  forms, 
several  other  from  a  root  locus,  and  the  last  one  or  two  by  direct  calculation. 

All  the  developments  of  this  section  are  carried  out  for  the  general  case 
of  an  Nth-order  closed- loop  pole,  but  at  a  convenient  step  the  results  for  the 
specialized  case  of  a  first-order  pole  are  often  given. 

To  furnish  a  concrete  example,  all  of  the  methods  developed  in  this  section 
have  been  applied  to  a  simple  third-order  system.  The  details  of  the  computa¬ 
tions  are  given  in  the  Appendix.  The  example  covers  the  use  of  the  techniques 
for  closed- loop  poles  which  are  real  and  first-order,  conplex  and  first-order, 
and  real  and  second-order.  The  results,  summarized  in  Article  E,  are  also  used 
to  illustrate  the  relative  accuracy  of  the  various  methods. 

A.  DIRECT  CALCULATION 


Of  the  many  methods  presented  in  this  report  for  the  evaluation  of  modal 
response  coefficients  and  gain  sensitivities,  only  the  two  given  in  this  article 
can  readily  yield  values  to  any  given  degree  of  precision.  All  the  other  tech¬ 
niques  involve  a  graphical  construction  or  are  based  on  certain  approximations. 
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Either  of  the  direct-calculation  isethods  could  be  added  to  a  digital  corrputer 
program  for  loop-closing  calculation  to  provide  the  modal  response  coefficients 
and  sensitivities. 


In  Section  II  the  modal  response  coefficients  and  gain  sensitivities  were 
given  by,  among  other  formulas. 
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Evaluation  of  Eq  125  is  particularly  convenient  when  the  numerator  and  denominator 
of  G  are  known  polynomials  in  s-  Expressing  G  as 
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From  Eq  125  and  125  it  is  obvious  that  for  N  =  1 


If  K  is  greater  than  one,  higher  derivatives  of  G  must  be  considered,  such 
as  a^G/as^  (Eq  1 26) .  Repeated  differentiation  of  Eq  1 26  combined  with  the 
requirement  that  the  first  (W  -  1 )  derivatives  of  G  with  respect  to  s  vanish  at 
s  =  -qjj^,  shows  that 
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For  H  =  1,  Eq  123  and  132  shov  that 
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which  is  identical  to  Eq  59. 
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For  N  greater  than  one  we  must  again  consider  higher  derivatives  of  G,  such 
as  Eq  133-  Repeated  differentiation  combined  with  the  vanishing  of  the  first 
(N  -  1)  derivatives  shows  that 
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The  coiroutation  of  Q^jj  or  by  Eq  1 37  can  easily  be  set  up  on  a  standardized 
work  sheet  and  would  probably  be  easier  to  use  in  a  digital  computer  program 
than  Sq  1 50-  It  is  also  siii5>le  to  show  that  Eq  1 57  is  valid  if  any  or  all  of 
the  terms  of  G  are  written  in  Bode  form- 


B-  ROOT-LOCUS  IffiTHODS 


Three  methods  of  using  a  root-locus  plot  to  estimate  the  modal  response 
coefficients  and  gain  sensitivities  are  described  in  this  article.  The  first 
two  are  approxinate  techniques,  while  the  accuracy  of  the  third  is  limited  only 
by  the  exactness  of  the  plot  and  the  required  measurements. 


The  first  method  requires  the  locations  of  the  closed-loop  pole  for  two 
slightly  different  gains,  such  as  the  segment  of  a  locus  shown  in  Fig.  5a* 
Using  finite  increments  as  approximations  to  differential  changes,  Eq  11 8  and 
122  give 
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Equation  138  can  also  be  utilized  if  the  change  in  the  closed- loop  pole  is 
obtained  by  some  other  technique,  such  as  from  |  plots  and  root  decomposition. 


The  second  perturbation  method  is  obtained  by  considering  /c  to  be  a  complex 
quantity.  The  normal  root  locus  is  then  a  graph  of  the  closed-loop  pole  loca¬ 
tions  for  K  real.  Now  consider  a  small  perturbation  in  the  phase  of  K.  The 
closed-loop  pole  must  then  be  perturbed  a  small  distance  normal  to  the  conven¬ 
tional  root  locus,  and  the  phase  perturbation  of  (G/k)  must  be  minus  that  of  k, 
see  Fig.  5b.  Consequently,  for  perturbations  normal  to  the  root  locus,  A/c/k  =  -jAS 
and  Eq  1  38  becomes 


S 


i 

K 


A© 


(159) 


where 


A0  =  phase  change,  in  rad,  of  G/k 


Note  that  for  N  greater  than  one  the  locus  is  not  a  single  curve  at  s  =  -qjj^,  but 
is  the  junction  of  N  incoming  and  N  outgoing  branches.  In  this  case  the  Aqj^  of 


Eq  159  should  be  picked  midway  between  any  incoming  branch  and  an  adjacent  out¬ 
going  branch. 

The  third  method,  which  will  be  called  the  vector  method,  has  the  advantage 
that  its  accioracy  is  limited  only  by  the  precision  of  the  graph  and  the  necessary 
measurements.  Of  the  graphical  techniques  described  in  this  section,  this  one 
should  generally  be  the  most  accurate.  Its  primary  disadvantage  is  that  it 
requires  a  coDiplete  set  of  conpatible  closed-loop  poles,  i.e.,  all  the  closed- 
loop  poles  for  a  particular  value  of  gain,  and  knowledge  of  the  value  of  that  gain. 


From  Eq  97  and  122, 


m+n 

(1  +  k6°)  n  (q.  -  q.) 


(l40) 


The  right  side  of  Eq  l4o  is  a  constant  times  a  function  which,  with  zeros  at  the 
open- loop  zeros  and  poles  at  all  the  closed-loop  poles  except  -q^,  is  evaluated 
at  s  =  -qj_.  Thus,  Eq  l4o  can  be  quickly  evaluated  using  the  normal  graphical 
root-locus  techniques  with  a  Spirule  or  similar  device.  The  principal  quantities 
for  a  sample  system  are  shown  in  Fig.  6. 

C.  METHODS  USING  OPEN-LOOP  BODE  AND  |  PLOTS 

The  derivative  ^%/Ss^,  and  consequently  the  modal  response  coefficient  and 
gain  sensitivity,  can  be  determined  from  open- loop  Bode  or  |  plots.  A  |  plot  is 
a  plot  of  the  amplitude  and  phase  of  G  as  a  fuuiction  of  p  for  constant  5,  where 


s 


{i  +  iVi  - 


(l4l) 


The  conventional  Bode  diagram  is  then  a  special  case  of  the  |  plot  (|  =  O) . 
Because  |  plots,  particularly  for  ^  =  ±1 ,  are  valuable  analytical  tools  for 


determining  the  closed-loop  poles,  it  is  highly  desirable  to  be  able  to  determine 
modal  response  coefficients  and  gain  sensitivities  from  the  |  plots. 


With  I  held  constant. 


3G  _  gSluG  _  g  d  In  u  3  In  G 
^  ~  5s  ~  ds  d  In  [1 


and  from  Eq.  1 4l , 


d  In  _  2 

ds  ~  s 


Equation  l42  can  now  be  written 


5G  ^  6  5  In  G 
5s  s  5  In  p 


The  second  derivative  is  then 


5s^ 


G_  5^  In  G  5  3n  G  5  ,Gn 
s2  5(ln  p)2  SIFp  5s 


Repeated  differentiation  combined  with  the  vanishing  of  the  first  (N 
tives  at  s  =  -q^  gives 


Representing  G  in  terms  of  its  eunplitude  and  phase,  i.e.. 


(142) 

(145) 

(144) 

(145) 

1 )  deriva- 

(1'^6) 


allows  us  to  rewrite  Eq  l46  as 


-1  #  In  A  ^ 

_^(lii  h{ln  }i)® 


-1  r  1  ^  log  A  ^  j 

L(2-3026)^"''  a(log  ixf  (2.3026)® 

_ ^1 _  [%  j%  1 

(-qi)®(2.3026)®-'  h  I5I.93J 

s — q. 


(l48) 


Therefore  for  N  =  1 , 


4  =  Qi  = 

and  in  general 


S 


i 

K 


-NI(2.5026)®"’'  (q^)® 

% 

20  151.95]  3..,,. 


(1‘^9) 


(150) 


Thus,  for  a  first-orier  pole  the  modal  response  coefficient  and  gain  sensi¬ 
tivity  can  be  determined  from  the  slopes  of  the  amplitude  and  phase  curves  of  a 
I  plot  at  s  =  -qj^.  If  q^  is  purely  imaginary,  a  conventional  Bode  diagram  is 
used;  if  qj^  is  real,  a  Siggy  =  ±1 )  diagram  is  used  and  the  q)  derivative  is 


zero.  For  a  complex  pole^  a  |  plot  or  the  shifted  Bode  diagram  method,  to  he 
described  shortly,  must  be  used. 


For  multiple-order  poles  the  higher  derivatives,  Ajj  and  (pjj,  must  be  evalu¬ 
ated.  This  might  be  done  by: 

1 .  Measuring  the  slopes,  A-j  and  cp^ ,  from  a  |  plot 

2.  Plotting  these  slopes  versus  log  p. 

5.  Measuring  A2  and  92 

4.  Plotting  A2  and  9^  versus  log  p 

5.  Repeating  the  process  until  Ajj  and  qjjj  are  found 

An  alternate  method  for  second-oider  poles  is  to  measure  the  radius  of  curvature 
of  the  I  plots.  The  radius  of  curvature  of  a  curve  defined  by  y  =  y(x)  at  a 
point  where  dy/dx  =  0  (Sg/Ss  is  zero  at  s  =  -q^)  is  given  by 


R 


(151) 


The  radius  of  curvature  can  be  used  to  determine  A2  and  92,  but  care  must  be 
exercised  to  properly  account  for  the  vertical  and  horizontal  scales  of  the  plot. 
A  simple  method  of  computing  A2  or  92  is  to: 

1  .  Measure  the  radius  of  curvature  (at  s  =  -q^^^)  in 
the  units  of  the  vertical  scale,  i.e.,  db  or  deg 

2.  Divide  the  vertical  scale  by  the  horizontal  scale, 
e.g.,  (db/in.)/(dec/in.)  or  (deg/in.)/(dec/in.) 

5.  Square  the  ratio  of  step  2  and  divide  by  the  radius 
of  curvature  measured  in  step  1 

An  example  of  both  methods  for  a  second-order  pole  is  given  in  the  Appendix. 

An  alternative  to  the  use  of  |  plots  for  complex  roots  is  to  use  shifted  Bode 
diagrams.  The  basic  idea  is  easily  understood  if  one  remembers  that  the  shape  of 
a  root  locus  is  independent  of  the  location  of  the  origin  of  the  coordinates. 
Consequently,  the  imaginary  axis  may  be  shifted  an  arbitrary  amount  without 
changing  the  locus  or  the  gain  sensitivity.  If  the  axis  were  translated  so  that 
it  passed  directly  through  a  complex  pair  of  closed-loop  poles,  the  gain  sensi¬ 
tivity  for  either  of  these  two  poles  could  be  obtained  from  a  Bode  diagram  and 


1.1 


Eq  150.  Since  Boie  templates  are  readily  available,  and  the  Bode  diagram  is  a 
well-hnoiin  analytical  tool,  the  shifted  Bode  diagram  may,  at  times,  be  easier 
to  use  than  the  §  plot. 

In  constructing  the  Bode  diagram  the  real  parts  of  all  the  poles  and  zeros 
oust  be  translated  by  the  proper  increment.  For  poles  or  zeros  on  the  real  axis 
this  involves  only  an  addition  or  subtraction.  For  complex  pairs  the  shift  will 
change  both  the  natural  frequency,  a>,  and  damping  ratio,  of  the  roots.  It  is 
easily  shown  that  if  the  iraginary  axis  is  shifted  to  the  left  (in  the  direction 
of  the  ne^tive  real  axis)  a  distance  d  =  Re(q^)  ,  the  natural  frequency,  cc^, 
aM  ratio,  ^3,  of  shifted  poles  or  zeros  are 


€0^  =  -  2^a)d  +  d^ 


r  = 

as 


!!<d  -  d 


fto  -  d 
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-  2^d  +  d^ 


For  equal  to  zero,  i-e.,  d  = 


=  (ay-i  - 


(155) 


A  third  method  using  |  plots  is  especially  appropriate  when  the  decomposition 
loop-closure  technique  is  applicable,  although  it  also  applies  in  general.  The 
basic  idea  is  to  find  two  sets  of  closed-loop  break  points  for  neighboring  open- 
loop  gains.  From  these  all  of  the  due  to  the  gain  change  can  be  estimated. 
Because  AC/K  =  Ac/ k,  Eq  1 58  can  be  rewritten 


= 


This  technique  is  the  Bode  diagram  version  of  the  gain  perturbation  technique 
discussed  in  connection  with  root-locus  methods.  Unfortunately,  the  method  is 
an  approximation,  i.e.,  differentials  are  replaced  by  increments,  and  may  have 
substantial  errors  if  appreciable  gain  changes  are  used  to  develope  the  basic 


D.  METHOD  USIHG  CLOSED-LOOP  BODE  ASYMPTOTES 


The  amplitude  asymptotes  of  the  closed- loop  Bode  can  frequently  he  readily 
obtained  hy  the  methods  of  Ref.  1 .  These  asyaiptotes  can  he  used  to  estimate  the 
modal  response  coefficients  and  gain  sensitivities.  This  is  only  an  approximate 
technique,  and  for  accurate  results  the  closed- loop  pole  must  he  widely  separated 
in  frequency  from  all  other  closed-loop  poles  (except  its  complex  conjvigate)  and 
from  all  the  open- loop  zeros. 


From  previous  developments  we  have 


n  (s  +  zjj) 

,  _  K  k=1 _ 

"  (1  +  /<50)  ,  V 

n  (s  +  q]j.) 


Qiw  = 


K  n  (zk  -  Gi) 

k=1 _ 

m+n 

(1  +  K&S)  n  (qjj.  -  q^) 
k=1 


If  -qj_  is  a  real  root,  the  products  of  Eq  156  can  he  approximated  hy 


(155) 


(156) 


n  (2k  -  <ii)  =  (-qi)^"  n '  (zk) 

k=1  k 

m+n  Xn  (m+n-W-Xg) 

n  (qk  -  qi)  =  (-qi)  n  (q,,) 

k=1  k 

•k^i 


(157) 


where 


n  (zk) 

k 


(m+n-K-Xg) 

n  (qk) 

k 


=  number  of  zeros  of  smaller  magnitude  than  q^^ 

=  product  of  the  (n-X^.)  zeros  of  larger  magnitude 
than  q^ 

=  number  of  closed- loop  poles  of  smaller  magnitude 
than  q^ 

=  product  of  the  (m+n-N-Xq)  closed-loop  poles  of 
larger  magnitude  than  q^^ 
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Combining  Eq  156  and  157, 


Q- 


=  = 


■lU 


X  “X«  h^Xt; 

K(-q.j)  ^  n  (Zfc) 

_  k  _ 

(m+n-N-Xq) 

(1  +  KbS)  n  (<3.k) 

k 


(158) 


From  Eq  155  it  can  be  seen  that  the  amplitude  asymptote,  of  the  closed 


loop  ^rc  Bode  diagram  at  a  frequency  equal  to  |qj^|  is 


asy 


1  +  /cbg 


X  n-X„ 

kil "  n  kki 


X„+N 


(m+n-K-Xa) 


(159) 


n 

k 


Note  that  several  terms  in  Eq  158  and  159  are  identical  except  for  the  absolute 
value  sign.  Consequently  these  equations  can  be  combined  to  give 


,  ,  .N  . 

(-1)  (<li)  A, 

(-1)  (“<1^)  ^asy 


asy 


q^  real  and  >  0 

(160) 

q^  real  and  <  0 


where  X^  =  number  of  closed-loop,  non-minimum-phase  zeros  and 

poles  of  larger  magnitude  than  q.- ,  plus  one  if 

k/(i  +  Keg)  <  0 


Note  that  because  complex  roots  always  occur  in  pairs,  the  powers  of  (-1)  in 
Eq  160,  (X2+Xq+Xu)  or  X^^,  can  be  considered  as  the  number  of  open-loop  zeros  and 
closed-loop  poles  which,  in  the  s-plane,  lie  on  the  real  axis  and  to  the  right 
of  -q^i^,  plus  one  if  k/(1  +  KSg)  <  0. 

If  -q^^  is  a  complex  pole,  Eq  158  must  be  modified  slightly  because  the 
complex  conjugate  pole  has  the  same  magnitude  as  q^.  Then  Eq  158  should  be 
written 


ims  IS  se±aoin  warrantea  oecause  tne  quicK-view  aavantag< 
and,  even  when  refined,  the  technique  cannot  conipete,  on 
the  others  discussed. 

E.  EXAMPLE 


To  illustrate  the  techniques  described  in  the  previo 
section,  all  of  the  methods  have  been  used  to  comnute  th' 


A  root-locus  plot  for  this  system  is  shown  in  Fig.  7.  The  calculations  were 
carried  out  for  two  different  values  of  gain,  K.  The  high-gain  case  was 
selected  to  illustrate  the  procedures  for  complex  closed- loop  poles,  and  the 
low-gain  case  illustrates  the  procedures  for  second-order  closed- loop  poles.  In 
addition,  the  actual  shifts  in  the  closed-loop  poles  due  to  the  gain  change  are 
conpared  with  the  values  predicted  by  both  the  high-gain  and  low-gain  sensitivities. 

For  the  high-gain  case,  the  gain  was  selected  to  produce  a  complex  pair  of 
closed- loop  poles  with  a  damping  ratio  of  -^272.  For  this  situation: 

K  =  2.070 

=  5-099 

qg  =  0.450(1  -  j) 

=  0.450(1  +  j) 

For  the  low-gain  case,  the  gain  was  selected  to  produce  a  second-order 
closed- loop  pole.  For  this  condition: 

K  =  1.128 

q^  =  5.055 

qg  =  0.472  (second-order  pole) 

The  values  of  gain  sensitivity  which  were  computed  by  e^'ch  method  are  listed 

in  Table  I.  The  values  of  for  the  high-gain  case  are  not  listed  because  they 

2 

are  the  complex  conjugate  of  S/<.  Wo  attempt  will  be  made  to  form  any  general 
conclusions  on  the  relative  accuracies  of  the  various  methods  except  to  note 
that  the  values  obtained  by  the  root-locus  vector  method,  which  would  be  expected 
to  be  the  most  accurate  graphical  technique  for  this  example,  agreed  very  well 
with  the  exact  values  obtained  by  the  direct-calculation  methods. 

To  get  some  indication  of  the  accuracy  of  predicting  the  changes  in  closed- 
loop  poles  from  the  gain  sensitivity,  we  can  compare  the  actual  shifts  between 
the  high-  and  low-gain  cases  with  the  changes  predicted  by  the  gain  sensitivities. 
This  comparison  is  summarized  in  Table  II.  Considering  that  the  ratio  of  the  two 
gains  is  more  than  1 .8,  the  agreement  is  generally  quite  good.  The  one  relatively 
poor  case  is  the  shift  in  qg  predicted  by  the  high-gain  sensitivity.  The  pre¬ 
dicted  change  in  the  imaginary  part  of  q2  is  only  half  of  the  actual  shift. 


62 


i^ 


a 

H  m 


a 

H  m 


o  I  a 
a  I  ei 


H  I  M 

cQ  I  a 


w  w 

a  H 

Eh  H 


o 

O 

KN 

00 

o 

VO 

trv 

ON 

ON 

IfN 

00 

ON 

VO 

5\ 

-d" 

KN 

KN 

• 

• 

• 

• 

• 

« 

« 

• 

• 

o 

o 

o 

o 

o 

o 

o 

o 

o 

•<T> 

TJ 

•o 

*o 

1 

1 

1 

t 

1 

i 

< 

1 

1 

t 

a 

o 

f— 

o 

m 

c^ 

t- 

KN 

If 

kS 

J- 

IfN 

a 

a 

o 

o 

O 

o 

O 

O 

d 

a 

d 

d 

d 

d 

d 

d 

1 

1 

1 

1 

1 

1 

■ 

i 

On 

ON 

CVJ 

KN 

VO 

KN 

ITN 

VO 

<3n 

ON 

ON 

03 

CVJ 

CJN 

ON 

o 

1^- 

O 

o 

O 

o 

O 

f— 

O 

rN 

« 

• 

• 

rs 

• 

o 

• 

o 

• 

O 

5h 

O 

•P 

(d 

S 

(U 

iH 

p 

CQ 

s 

CQ 

d 

0) 

3 

CQ 

a 

0) 

p 

■P 

S 

s 

a 

o 

ft 

Oj 

p 

0 

dj 

o 

p 

o 

> 

O 

•P 

■d 

<u 

•H 

p 

rH 

p 

p 

ft 

p 

p 

cd 

CQ 

d 

p 

s 

Ti 

cd 

P 

o 

cd 

d 

'h 

P 

P 

0) 

<u 

P 

CQ 

oJ  m 

o 

5h 

d 

> 

ft 

•d 

P 

Cd 

<u 

d 

p 

o 

o 

d 

u  > 

a 

p 

p 

CQ 

m 

p 

Q> 

O  -H 

o 

fH 

0) 

CQ 

CQ 

p 

P  P 

p 

<u 

Pi 

<D 

d 

•d 

0) 

(d  d 

p 

Pi 

p 

O 

p 

<u 

ft 

■p 

Jh  > 

<U 

o 

O 

•d 

p 

'H 

(U  P 

M 

P 

m 

p 

(d 

ft 

p 

tH 

§  ^ 

n 

P 

% 

o 

CQ 

p 

p 

p 

ft 

3  (1) 

d 

Cd 

0) 

a 

cd 

& 

a  -d 

03 

o 

a 

> 

UP 

u/> 

03 

O 

c 

•d 

o 

P 

•H 

cd 

0) 

p 

(1) 

>d 

o 

m 

d 

o 

o 

W 

ft 

P 

CQ 

O 

cd 

d 

ft 

o 

o 

CJ 

o 

p 

o 

o 

(0 

1 

p 

p 

P 

p 

'd 

o 

1 

o 

d 

0) 

p 

p  P 

(0 

!h 

o 

(U 

ft 

o 

•H 

o 

ft 

P 

ft 

a 

o 

UP 

o 

H 

M  « 


^  K^  lA 


Predicted  by  high-gain  sensitivity 
Predicted  by  low-gain  sensitivity 


SECTION  IV 


MULTILOOP  SYSTEMS 

A.  MODAL  RESPONSE  COEFFICIENTS 

The  modal  response  coefficients  for  a  multiloop  system  can  he  derived  in  a 
manner  quite  similar  to  that  used,  in  Section  II  for  single-loop  systems.  The 
familiar  egression 


also  applies  to  multiloop  systems,  hut  Gj.g  is  now  a  more  complicated  function. 
In  general,  Gj.q  can  he  written  as 

F 

^rc  1  +  Fo 


where  F-j  and  F2  are  sums  of  terms  which  contain  the  transfer  functions  of  the 
various  individual  elements  in  the  system. 


As  in  Secoion  II,  define 


qf  (s)  =  (s  +  qi)\ 


(s  t  qi)V 


1  +  Fc 


(1  +  F2)Q^jj(s)  =  (s  +  q,^)\ 


Repeated  differentiation  of  Eq  167  shows  that 


S^Fp 

Ss^  ~ 

•-  -'s=-a-- 


1  <  k  <  N  -  1 
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A  relationship  for  the  summation  for  the  modal  response  coefficients  can 
also  he  derived  in  a  manner  analogous  to  that  of  Section  II.  If  there  are  no 
multiple-order  closed- loop  poles,  and  if  the  number  of  closed-loop  poles  is 
greater  than  the  number  of  closed- loop  zeros,  we  can  write 


'rc 


,  E 

s  +  a,- 


(170) 


-1 


Multiplying  both  sides  by  s,  and  letting  s  approach  infinity,  gives 


EOi 


if  niamber  of  closed- loop  poles  =  1  +  number 
of  closed- loop  zeros 

0  if  number  of  closed-loop  poles  >  2  +  number 
of  closed- loop  zeros 


(171) 


where  Kqj^  is  the  closed- loop  root-locus  gain  of  the  system.  This  is  directly 
analogous  to  the  result  obtained  for  the  single-loop  case. 


When  there  are  Nth-order  closed-loop  poles,  a  result  analogous  to  that  for 
single -loop  systems  is 


E^i  E*^i-| 

i  i 


if  number  of  closed-loop  poles  =  1  +  number 
of  closed-loop  zeros 

0  if  number  of  closed- loop  poles  >  2  +  number 
of  closed-loop  zeros 


(172) 


B.  SENSITIVITIES 


For  multiloop  systems  the  closed- loop  poles  are  functions  of  the  gains,  as 
well  as  the  pole  and  zero  locations,  of  a  number  of  transfer  functions.  It  is 
then  desirable  to  know  the  sensitivity  of  a  closed- loop  pole  to  variations  in 
the  gain,  poles,  and  zeros  of  any  of  the  transfer  functions.  To  derive  these 
sensitivities,  consider  the  characteristic  equation  of  a  multiloop  system. 


s=-q. 


0 


u  y:); 


where  F2  is  a  summation  of  terms,  each  of  wnich  is  the  product  of  various  trans¬ 
fer  functions.  To  keep  the  discussion  completely  general,  the  transfer  function 
of  any  individual  element  in  the  system  will  he  denoted  by  Gj^.  Katurally  the 
functional  relationship  between  F2  (or  F^ )  and  any  individual  transfer  function 
depends  on  the  manner  in  which  the  various  elements  are  connected.  This  func¬ 
tional  relationship  can  be  determined  from  a  block  diagram  or  signal  flow  diagram 
of  the  system. 

Extending  the  notation  to  the  gains,  poles,  and  zeros,  is  the  root-locus 
gain  of  the  kth  transfer  function,  and  (-zj^.)  and  (-pj^.)  are  the  jth  zero  and 
pole  of  the  kth  transfer  function.  The  total  differential  of  F2  can  be  written 
(analogous  to  Eq  II6  and  11 7) 


where  nj^  =  number  of  zeros  of  kth  transfer  function,  Gjj 

mk+n]j^  =  number  of  poles  of  kth  transfer  function 
L  =  number  of  transfer  functions 


Note  that  as  per  Eq  1 68,  the  first  (N  -  1 )  derivatives  of  F2  with  respect  to  s 
have  been  set  equal  to  zero. 
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As  in  the  single-loop  case,  we  will  define  the  sensitivities  so  that 


Fl  /  .  a(c.  .  V“k  .  0’^“ 


(176) 


Equating  like  coefficients  of  Eq  175  an<i^  17^  gives 


L  Js=-qi 


s=-qi 


(177) 


(178) 


^k 

— fj- 

§G^ 


(179) 


^.1 

71S 

k  ^ 


(180) 


If  the  transfer  functions  are  written  in  root-locus,  Bode,  or  adxed  form. 


(181) 
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Consequently, 


(182) 


s=-q. 


6^2 


(185) 


|s=-q^ 


For  zeros  or  poles  which  are  written  in  root-locus  form. 


%  - 


(184) 


(185) 


Likewise,  for  zeros  or  poles  which  are  written  in  Bode  form. 


^i 

zkj(zkj  -  <iiy 


(186) 


PkjWi  -  Pk.> 


(187) 


Two  important  features  of  the  above  sensitivity  ratios  are  worth  mentioning. 


First,  because  the  terms  of  F2  will  normally  be  products  of  transfer  fionctions 


r.  -• 
r: 

-  •  .  *  • 


«  •  •  •  »  •  • 


-  - 


^  *’  s'  •.*  v*  s"  •. 


to  the  first  power,  Gjj;(3F2/3G]j.)  will  simply  he  the  sum  of  all  the  terms  of  F2 
which  contain  Gj^.  Thus,  the  ratio  of  the  gain  sensitivities  for  k-  and  will 
be  the  ratio  of  the  terms  of  F2  in  which  Gj  and  Gj^^  appear,  evaluated  at  s  = 
Second,  the  ratios  of  zero  or  pole  sensitivities  to  the  gain  sensitivity  of  the 
same  transfer  function  are  identical  to  those  for  the  single-loop  case.  Conse¬ 
quently,  the  previously  derived  (Eq  78  through  8l )  sensitivity  relationships  in 
terms  of  the  damping  ratio  and  frequency  or  the  real  and  imaginary  parts  of  an 
open- loop  root  also  apply  to  the  multiloop  case. 

Because  the  normal  procedure  for  evaluating  the  closed-loop  poles  of  a  multi¬ 
loop  system  is  by  means  of  a  series  of  loop  closures,  the  above  equations  can  be 
utilized  to  form  a  relatively  simple  procedure  for  evaluating  the  sensitivities: 

1 .  Calculate  one  of  the  gain  sensitivities  by  applying  one 
of  the  methods  of  Section  III  to  the  final  loop  closure 

2.  Compute  the  other  gain  sensitivities  from  Eq  183 

3*  Conpute  zero  and  pole  sensitivities  from  Eq  78,  79,  80. 

81,  184,  185,  186,  or  187 

In  many  cases  step  2  of  the  above  will  be  the  greatest  source  of  computational 
difficulty.  Wo  simple  standard  procedures  can  be  established  for  this  operation, 
but  there  are  several  ideas  which  may  be  helpful.  One  of  these  is  to  make  use  of 
the  fact  that  at  s  =  -q^,  Fg  =  -1.  For  example,  if  F2  were 


F2  =  G1G2  +  GiGjG]|^  +  G']G2G^  + 


then 


( _ ^-1^2  \  ^  _\[ 

yGiG2  +  G-iG^Gij^  +  GiG2Gji  | 

's=-qi 


^162(1  +  Gj) 


1  +  G^Gj^  j 


s=-q. 


/G-jGjGij^  +  G-]  G2G^  +  G^G]j^^ 

/1  +G^G2\ 

^G^G2  +  G^GjGj^  +  G,G2G5j 

V 

s=-q. 


s  =  -q_. 


qi 


/  G^G3G^-^G^G^  \  R 

\^G^G2  +  G^GjGi^  +  G^G2gJ  " 

's=-q.  ^ 


+G^)' 


1  + 


s=-q_. 


••'-.’J 


m 


^ A 


A  Aj 


'  .  •  A 

-'.'•"a 

U-AiliuS 
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The  sensitivity  ratios  could  be  computed  directly  by  setting  s  =  in  the 
proper  equations,  but  for  coii5)lex  poles  graphical  methods  may  be  simpler.  The 
products  of  transfer  functions  could  be  evaluated  by  constructing  I  plots  or 
root- locus -like  plots.  From  the  known  zeros  and  poles  of  the  transfer  functions, 
a  I  plot  could  be  constructed  for  |  equal  to  minus  the  damping  latio  of  q^^.  The 
constixiction  would  only  have  to  be  accurate  for  the  frequency  equal  to  |  q^  | . 

The  an5)litude  and  phase  of  the  product  of  transfer  functions  can  also  be  evalu¬ 
ated  by  plotting  zeros  and  poles  on  an  s-plane  graph,  and  measuring  the  amplitude 
and  phase  at  s  =  -q^.  Once  the  products  have  been  evaluated,  the  summations 
with  unity,  which  are  sometimes  necessary,  could  be  done  nianerically,  from  a 
Nichols  chart,  or  graphically  (preferably  on  polar  graph  paper) . 


As  in  the  single-loop  case,  a  useful  check  on  the  sensitivity  calculations 
is  that  the  sum  of  all  the  zero  and  pole  sensitivities  is  one  if  the  transfer 
fimctions  are  entirely  in  root-locus  form  and  if  N  equals  one  (-q^^  is  a  first- 
order  pole).  This  can  be  shown  from  Eq  l82,  and  185. 


L  /"k  .  “k+“k  .  \  r  .  L  8Fp/°k  . 

iS  (j§  °  iS  ^(j?i  "kj  - 


jS  -  p  J 


The  derivative  8F2/8S  can  be  written  as 


(188) 


J  s=-q. 


\ 


Consequently, 


+  E  sj 


1 


(190) 


Sis  zeTationship  betireen  "tbs  gain  ssajsitivily  and  the  coaal  response 
coefficient  can  be  seen  directly  fjraa  Ba  169  and  l82. 


=  (-ir 


O^!) 


s=-qi 


For  a  snaltiple  iiroat-ontont  system  there  will  be  a  set  of  nodal  response  coeffi¬ 
cients  for  each  is^ut-output  cac^ination.  All  the  coefficients  can  be  e-valuated 
by  means  of  Sn  191  if  the  prt^r  msKrator,  F|,  for  each  innut-ontpat  pair  is  used. 

C.  EZffiSSIE 

altitnie  control  system  of  Bef.  2  will  be  used  as  an  exarrnle.  In  that 
^sfcea,  pitch  angle,  0,  and  altitude,  h,  are  fed  hack  to  the  elevator.  Fig-  8. 


Figure  8-  tJultiloop  Exasple 


In  Fig-  8,  Sg  and  Hg  are  the  airfrasoe  pitch  angle  and  altitude  transfer  functions, 
Yg  and  are  the  equalization  and  sensor  dynandcs,  and  Yg  is  the  elevator  servo 
dynamics.  From  Fig.  8  it  can  be  seen  that  the  closed-loop  transfer  function  is 


f.'-r  ■ 


Yh%Y5 

1  -f  YgHgYg  +  YgegYg 


O92) 


For  convenience,  the  notation 


G2  =  Yg  (195) 

^3  =  ^,05 


will  be  used.  The  form  of  the  transfer  functions  and  numerical  values  which  were 
used  are: 


G,  = 


0-2  = 


"  %)t '  s) 

(s  +  a^)s(s2  +  2^pa^s  +  a|)(s2  +  2^sp‘%pS  +  t>|p) 
% 


G3  =  K^Ag 


/  .  1  \ 

(""“Jr  ^tt) 

\  /  \  ‘'1/ 

The  characteristic  equation  of  this  system  is 


(l  +  G-jG2  +  G^G'^) 


s=-qi 


(195) 


The  technique  for  determining  the  closed-loop  poles  which  was  used  in  Ref.  2  was 
to  rewrite  the  characteristic  equation  as 


s--q^ 


Equation  196  can  also  be  written  as 


6  ^  ^  4)('  ■  '^3)1 


(196) 


the  complex  conjugates  of  the  computed  sensitivities 


"These  values  for  the  closed- loop  poles  are  slig 
given  in  Ref-  2.  The  values  given  above  were  obtain 
annroximation  to  the  closure  thain  was  done  in  Ref.  2 


The  firot  step  is  to  determine  one  of  the  gain  sensitivities.  The  simplest 
one  to  do  is  the  one  for  K^^  (the  gain  of  G-])  using  the  closure  defined  hy  Eq.  197* 
Any  one  of  the  techniques  of  Section  III  could  he  used  to  evaluate  ,  hut 
because  all  of  the  other  sensitivities  will  he  computed  from  ,  we  will  use 
the  exact  analytical  expression 


4,  ■  [f 


k  ^k-% 


5? 


(199) 


where  -pj^  and  are  the  poles  and  zeros  of  G^ 62/(1  +  62^3)  ^  ^*1  19^  19’^* 

This  computation  gives 


=  -(0.175  +  jo.492) 


=  -0.521  exp 


(‘Sj) 


(200) 


For  the  sake  of  brevity,  we  will  adopt  the  notation 


Aej^P  H  A  4 


(201) 


and  express  as 


Sj;  =  0.521  4  250.4  deg 


(202) 


Note  that  because  =  G-|  G2  =  G-j  (^F2/SGi  )  ^  is  equal  to  the  modal  response 
coefficient  for  h/h^- 

As  the  next  step,  and  S^  will  he  computed.  From  the  sensitivity  ratios 
of  Eq  183,  ^ 


SL  \  ^1^2  A 


(203) 


si  1^1  ^2/  \^l/ 

\  's=-q^  '  's=-q^ 


(204) 
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By  plotting  the  poles  and  zeros  of  &|  and  G2  on  an  s-plane  plot,  (tf]  G2/k-i  K2) 
can  be  measured  with  a  Spirule  or  similar  device.  Multiplying  by  the  gains  gives 


-T-  =  0.160  4  53-3  deg 

% 

Therefore, 

^K2  ""  (0.521)  (0.1 60)  4  (250.4  +  55.3)  deg 
=  0.0835  4  303.7  deg  =  0.0464  -  jO.0694 

Similariy,  by  plotting  the  poles  and  zeros  of  G5  and  1  /G^  , 


(205) 


(206) 


%  (h] 

~  V^l/ 


=  0.965  4  172.4  deg 


s=-q.- 


=  0.504  4  62.8  deg  =  0.230  +  jO.446  (207) 


From  Eq  203  and  204  it  can  be  seen  that  +  S j  should  be  zero.  Using 

the  coniputed  values  gives 


%  +  =  0.009  +  jO.023 

Considering  the  graphical  techniques  used,  the  agreement  is  quite  good. 

Now  the  pole  and  zero  sensitivities  can  be  computed  from  the  sin^jle  sensi¬ 
tivity  ratios  of  Eq  l84  and  185.  The  results  are  summarized  in  Table  III.  It 
was  noted  earlier  in  this  section  that  the  sum  of  all  the  pole  and  zero  sensi¬ 
tivities  of  each  closed- loop  pole  should  equal  one.  The  actual  sum  of  the 
sensitivities  of  Table  III  is  1 .009  +  j0.029.  This  agreement  is  excellent 
considering  the  inaccuracy  of  the  closed- loop  pole,  the  graphical  techniques 

used  in  the  calculations,  and  the  fact  that  the  numerical  computation  was  all 
done  by  slide  rule. 


K 


TRANSPER 

FUNCTION 


TABLE  III 

2ERO  AND  POLE  SENSITIVITIES  FOR  MULTILOOP  EXAMPLE 

SENSITIVITY 


0 

0.568  2^ 

-57.9  deg 

= 

0.448 

- 

jO.548 

15 

0.0553  -2^ 

67.0 

= 

0.0l4 

+ 

jo.055 

0.004-5  +  j0.0628 

0.535  A 

-36.4 

= 

0.431 

- 

jo.317 

0.0045  -  j0.0628 

0.609  A 

-38.8 

= 

0.475 

- 

00.582 

2.10  +  j3-7l 

0.1058  4 

2.0 

= 

0.106 

+ 

j0.004 

2.10  -  j3.7l 

0.1547  4 

127.9 

= 

-0.095 

+ 

j0.l22 

(0.0064) 

0.571  4 

i42.5 

= 

-0.452 

+ 

jO.347 

(19.2) 

0.0275  4 

-112.2 

= 

-0.010 

- 

jO.025 

(-19.2) 

0.0267  4 

73.0 

= 

0.008 

+ 

OO.O26 

35.0  +  035.7 

0.001654 

77.2 

= 

0.0004 

+ 

j0.00l6 

55.0  -  j35.7 

0.001704 

168.8 

= 

-0.0017 

+ 

jO.0003 

0.0045  +  j0.0628 

0.516  4 

136.0 

= 

-0.371 

+ 

CO 

• 

0 

•0 

0.0045  -  j0.0628 

0.590  4 

155.6 

= 

-0.4o6 

+ 

jO.426 

2.10  +  J5.71 

0.1025  4 

174.4 

= 

-0.102 

+ 

jO.OlO 

2.10  -  j3.7l 

0.1492  4 

500.5 

= 

0.075 

- 

jO.129 

(2.4) 

0.221  4 

4o.4 

= 

0.168 

+ 

j0.i43 

(0.0098) 

0.552  4 

515.1 

= 

0.591 

- 

jO.389 

(1 . 572) 

0.361  4 

24.1 

= 

0.530 

+ 

j0.i48 

•<*-  *7.  •  A  ^  ^  -  •L M  ^  VTrf  *1  _  •r^l ^  k  ^  iT-.^ ».  ■ 


Note  that  the  phugoid  (a^)  and  short-period  (ojsp)  poles  appear  in  both  Gi  and  G5 
Consequently,  the  total  sensitivities  of  these  poles  is  the  sum  of  their  sensi¬ 
tivities  when  considered  as  poles  of  Gq  and  of  G^.  It  is  also  desirable  to 
convert  the  sensitivities  for  the  complex-pair  poles  to  damping  ratio  and 
frequency  sensitivities  per  Eq  78.  The  final  sensitivities  of  all  the  variable 
parameters  are  summarized  in  Thble  IV.  The  gain  sensitivities  are  separated 
from  the  pole  and  zero  sensitivities  as  a  reminder  that  gain  sensitivities  are 
based  upon  fractional  variations  rather  than  absolute  variations. 

TABLE  IV 


MULTILOOP  EXAMPLE  SUMMARY 
q  =  0.287  -  jO.869 


GAIN  SENSITIVITIES 

Parameter 

Description 

Sensitivity 

Gain  of  altitude  feedback 

-0.175  -  jO.492 

Xe 

Gain  of  pitch-angle  feedback 

0.230  +  jO.446 

'^m 

Gain  of  elevator  servo 

0.046  -  jO.069 

■^h 

Gain  of  Hg 

-0.175  -  jO.492 

Gain  of  0g 

0.230  +  jO.446 

POLE  AND  ZERO  SENSITIVITIES 

Parameter 

Description 

Nominal 

Value 

Sensitivity 

Altimeter  lag 

15.0 

0.014  +  jO.055 

Lead  in  pitch-angle  feedback 

2.k 

0.168  +  j0.l45 

^m 

Damping  ratio  of  elevator  servo 

0.7 

-0.0019  -  j0.0024 

% 

Frequency  of  elevator  servo 

50.0 

-0.0018  +  j0.0028 

1/Thi 

First  zero  of  Hg 

0.0064 

-0.452  f  jO.347 

OJ 

Second  zero  of  Kg 

19.2 

-0.010  -  jo.025 

1/Th3 

Third  (unstable)  zero  of  Hg 

19.2 

-0.008  -  jO.026 

1/T01 

First  zero  of  0g 

0.0098 

0.391  -  jO.389 

l/Tsp 

Second  zero  of  65 

1 .572 

0.550  +  j0.l48 

Damping  ratio  of  phugoid  mode 

0.0714 

0.0084  +  jO.0033 

Frequency  of  phugoid  mode 

0.063 

0.0121  -  j0.0045 

^sp 

Damping  ratio  of  short-period  mode 

0.495 

-0.018  -  j0.028 

^sp 

Frequency  of  short-period  mode 

4.27 

-0.025  +  jO.025 

An  examination  of  Table  IV  reveals  a  number  of  interesting  bits  of  informa¬ 
tion: 

1  .  The  gain  sensitivity  of  and  Aj^  is  nearly  equal  to  the 
negative  of  that  of  and  Aq 

2.  The  gain  sensitivity  of  is  much  less  than  that  of  the 
other  gains 

5.  On  the  basis  of  the  same  percentage  change  in  parameters, 
the  dominant  mode  is  most  sensitive  to  o^,  03^,  1  /Tg  , 
l/Thg,  and  1/Th^ 

We  now  know  which  parameters  are  the  most  important  to  the  dominant  mode  of  the 
system,  and  can  estimate  the  variations  in  the  dominant  roots  which  could  be 
achieved  by  changing  the  gains.  By  combining  the  above  results  with  "approxi¬ 
mate  transfer  fimctions"  or  an  equivalent  technique,  which  can  relate  the  zeros 
and  poles  of  the  aircraft  transfer  function  to  the  stability  derivatives,  we  can 
estimate  the  changes  in  the  dominant  roots  due  to  variations  in  the  flight  con¬ 
ditions  and  due  to  uncertainties  in  the  inertial  and  aerodynamic  characteristics 
of  the  vehicle. 


SECTION  V 


SUMMARY 

A  good  many  derivationii  and  relationships  are  scattered  throughout  this 
report.  To  facilitate  the  use  of  the  material  contained  herein,  the  important 
results  have  been  summarized  in  a  number  of  tables  which  are  contained  in  this 
section. 

Table  V  lists  the  relationships  between  the  modal  response  coefficients  and 
the  system  responses  for  various  inputs.  Part  A  contains  the  relationships  for 
N  =  1 ,  and  Part  B  contains  the  general  relationships . 

Table  VI  lists  a  number  of  useful  identities.  Part  A  lists-  the  identities 
for  N  =  1 ,  and  Part  B  lists  the  general  identities . 

Table  VII  contains  the  sensitivity  ratios.  Part  A  is  the  ratios  for  zeros 
or  poles  which  are  in  root-locus  form.  Part  B  is  the  ratios  for  zeros  or  poles 
in  Bode  form. 

Table  VIII  summarizes  all  the  methods  for  computing  gain  sensitivities  which 
were  developed  in  Section  III.  The  table  lists  the  pertinent  equations  along 
with  some  of  the  advantages  and  disadvantages  of  each  method. 


TABLE  VI-A 


EDEMTITIES 


Single  Loop  vith  Unit  Feedback 


Qi  =  si 


m+n 

S  % 

i=1 


K  ,  m  =  1 
0  ,  m  >  2 


/c(bi  -  a-i  +  k)  ,  m  -  1  » 

m+n  Y  — 

D  Qi(li  =  -/c  ,  m  =  2  iti  ‘ii 

i=1 

0  ,  m  >  5 


1  ,  k  >  1 


1  +  K 


,  k  =  0  >  m  >  1 


0  ,  k  <  -1 


“  i=1  % 


m+n  Q. 
1=1 


53  S^.  +  ]^  Sp.  =  1  ,  all  zeros  and  poles  in  root-locus  form 


Multiloop 


Qi  = 


E  Qi  = 


number  of  closed-loop  poles  =  1  +  number  of 
CL  ’  closed- loop  zeros 

^  number  of  closed- loop  poles  >  2  +  number  of 


0  , 


closed-loop  zeros 


L  /  “k  “k+"k 

E  Z  ^  z  sj 

k=1  \j=1  j=1  H 


=  1  ,  all  zeros  and  poles  in  root- 

locus  form 


TABLE  VI-B 


IDEHTITIES 

GENERAL 


Single  Loop  with  Unit  Feedback 


S 


i 

K 


E  Qi 


K  ,  m  =  1 
0  ,  m  >  2 


Multiloop 


% 


=  (-1) 


N-1  r  1  qi 

Gi,(dF2/dGk) 

i-  ^i 


E  Qii 


^  number  of  closed- loop  poles  =  1  +  number  of 

CL  ’  closed-loop  zeros 

Q  number  of  closed-loop  poles  >  2  +  number  of 
^  closed-loop  zeros 


TABLE  VII-A 


SENSITIVITY  RATIOS 
ZEROS  OR  POLES  IN  ROOT-LOCUS  FORM 


"  -  % 

sj. 

^^k 

- 

^03 

*±2(a3-Cqi)Si 
qf  -  2^axij^  + 

S^ 

*+2axi^S^ 
q?  -  2^0X1^  +  0)^ 

4 

*±2(a  -  qi)S^ 
q?  -  2aqj|^  +  a^  +  b^ 

Sb 

*±2bS^ 

q?  -  2aq^  +  a^  +  b^ 

Use  upper  sign  for  zeros,  lower  sign  for  poles. 
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TABLE  VII-B 


METHODS  OF  COMPUTIKC 


Nuirierator  s.'s!  der.sr.ir^tor 
derivatives 


Open-loop  Bode  and  Successive  slopes 
i  plots 


Radius  of  curvature 


Gain  perturtatlo 


Closed-loop  Bode  I  Anplllude  esjarptoies 


L  OS  osj 


«  n  -  q?) 
.1=1 


B+n 

(1  +  k6§)  n  (qj  -  qj) 
j/i 


03s  •  Vo?  -  S^cul  +  d^ 


•  a>. 


(-1)  “  *  ' 


,  real  ar.d  >0 

,  real  ar.d  <0 


,  '+>'u/‘'i*  \^s‘^q  *osv 


S*  -  (-i)'‘''’k 


(-iA: 


/  I  >* 


I  vertical  sc 
1  horizontal  < 


TABLE  VIII 


I-ETHODS  OF  C014PUTING  GAIK  SENSITIVITIES 
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qj  cor.plex 


REyARKS 


J .  S<=plest  nu=erlcol  seUiod  for  loa-onJer  systess- 
2.  Requires  open-loo?  nuseretor  and  densalnator  polynoalals- 


1.  Sisplesi  nuaericnl  aelhod  for  hlgh-order  systeas. 


]•  Siaplest  root- locus  ecthod.  tJ)  easily  obtained  vlth  Splrule. 


1 .  OrdirArily  the  aost  accurate  graphical  technique. 

2.  Requires  coaplete  set  of  corpatible  closed-loop  poles. 


1 .  Heed  only  the  portion  of  the  Bode  or  {  plots  about  |c|  •  I..  ■ 

2.  Accurate  acasureaents  of  slopes  or  radii  of  curtature  are 
frequently  difficult  to  obtain- 


).  Method  Is  siaple  and  speedy  when  decoaposition  technique  is 
applicable. 

2.  Sensitivities  for  all  closed-loop  poles  can  be  coaputcd  froa 
one  gain  change  when  decoaposition  technique  Is  applicable. 


1.  Method  is  oidinarily  siaple  and  provides  quick  overview  when 
closod-loop  Bode  asysiptotie  plot  is  available. 

2.  T.hc  cloocd-ioop  pole  must  be  widely  separated  in  frequency 
from  all  zeros  (except  dipole  zero,  -zj,)  and  f.-on  all  other 
Closed-loop  poles  except  its  coaplex  conjugate;  to  correct 
for  dipole,  multiply  sensitivity  by 

('  -  <  ^5ll  ^  I’ll  <  l^hl 
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APPENDIX 
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SINGLE-LOOP  EXAMPLE 

The  vise  of  all  the  methods  for  con^iuting  gain  sensitivities  and  modal  response 
coefficients  will  he  illustrated  for  a  simple  exanple.  The  sample  system  has  an 
open- loop  transfer  function  given  by 

^  =  s(s  +  lHs  +5) 

The  computations  will  be  given  for  two  different  gains.  The  high-gain  case 
illustrates  the  procedure  for  complex  closed- loop  poles,  and  the  low-gain  case 
is  an  example  of  the  computations  for  a  second-order  closed- loop  pole. 


A.  HIGH  GAIN 

For  this  case  the  gain  is  set  so  that  a  pair  of  complex  closed- loop  poles 

with  damping  ratio  of  ■^^/2  exists.  For  this  situation  the  important  parameters 

K  =  31 -  156  =  2.070 

=  ^25  =  5.099 

q2  =  5  (6  -  y/25)(l  -  j)  =  0.450(1  -  j) 

qL5  =  ^  (6  -  y/^)0  +  j)  =  0.450(1  +  j) 

1 .  Direct  Calculation 


Using  Eq  127  for  the  numerator  and  denominator  derivatives  method, 


% 

From  Eq  A-1 ,  a 

3a 

P 

SP 


0 

•z  o 

s-^  +  6s  +  5s 
3s^  +  12s  +  5 


(A-2) 


2.  Root-Locus  Methods 


The  root  locus  of  the  exainple  system  was  shown  in  Fig.  7*  To  estimate  the 
gain  sensitivity  from  Eq.  138,  a  position  on  the  locus  of  -0.4l  +  jO.76  was 
chosen.  The  measured  gain  at  that  point  was  5.89.  This  gives  an  estimate  of 


o2  .  '^2  ^  [O.^H  -  jO.76  +  (-0.45  +  jO.45)] 

K  -  -  '^•07  5.89  -  2.07 

=  -0.045  -  jO. 355  =  0.3564  262.7  deg 

which  has  an  a]i5)litude  error  of  28  percent  and  an  angular  error  of  2  deg. 
For  the  point  -5.42  was  selected,  and  the  gain  was  IO.O6.  Then 


o1  ^  2.07(5.42  -  5.10) 
K  10/06  -  2.07 


0.085 


Despite  the  fact  that  the  gain  was  increased  to  nearly  5  times  its  original  value, 
the  estimate  is  within  1  3  percent  of  the  exact  value . 


To  obtain  estimates  from  Eq  1 39;  perturbations  normal  to  the  locus  were  con¬ 
sidered.  For  the  point  -0.64  +  o0.47  was  selected,  and  the  measured  phase 
change  was  -24  deg.  Then 

j[0.64  -  jO.47  +  (-0.45  +  jO.45)] 

-24/57.5 

-0.0477  -  jO.454  =  0.457  4:264  deg 

This  estimate  is  fairly  good;  the  amplitude  error  is  7  percent  and  the  angular 
error  is  less  than  1  deg. 


St  = 


A0 


For  q^  the  point  -5. 10  +  jO.10  was  chosen.  The  phase  change  was  47.5  deg, 
so 


S 


1 

K 


j(j0-10)(57.5) 


0.121 


This  estimate  is  in  error  by  27  percent. 


These  values  are  extremely  close  to  the  exact  values.  Both  amplitude  errors  are 

2 

less  than  1  percent,  and  the  angular  error  of  is  less  than  1  deg. 

5.  Methods  Using  Open-Loop  Bode  and  |  Plots 

A  I  plot  of  the  example  for  |  =  -O.7  is  shown  in  Fig.  A-1  .  A  small  error 
has  been  introduced  into  the  solution  from  the  use  of  the  |  =  -O.7  teirplates 
because  the  actual  value  of  |  at  q2  “  -0.707*  The  slopes  which  were 

measured  from  this  figure  are 

A^  =  -17  db/dec 


-129  deg/dec 


Then  from  Eq  l49 


[-0.450  +  JO.450] 


20  "  ^ 


.  129 
3 


=  -[0.0541  +  00.490]  =  0.491  4  266  deg 


This  is  within  1  percent  of  the  correct  amplitude  and  2  deg  of  the  correct 
angle.  The  extremely  small  error  in  the  aii?)litude  must  be  considered  fortuitous 
as  the  errors  in  measuring  the  slopes  will  normally  be  higher  than  this. 


For  we  must  use  a  ^  =  -1  or  Siggy  plot,  see  Fig.  A-1 .  The  measured  slope 
from  this  plot  is  -83O  db/dec,  but  the  root  is  too  close  to  the  open- loop  pole 
to  get  an  accurate  value.  A  better  estimate  can  be  obtained  by  noting  that  the 
two  low  frequency  roots  contribute  about  -4o  db/dec,  and  the  contribution  of  the 
pole  at  -5  can  be  approximated  by  (Ref.  1 ,  Table  III-A) 


.20^1^ 

:).io  -  5 


=  -1020  db/dec 


Therefore, 


A^  =  -1060  db/dec 


20(-5.10) 

=  '  -1060  ■  = 


This  is  within  2  percent  of  the  exact  value. 


Instead  of  using  the  i  plot  for  S^,  we  might  use  the  shifted  Bode  plot. 
Shifting  the  imaginary  axis  0.450  to  the  left  so  that  and  q^  are  on  the  axis 
puts  the  open-loop  poles  at  0.450,  -0.550^  and  -4.550*  The  shifted  diagram  is 
shown  in  Fig.  A-2.  The  measured  slopes  are 


A-]  =  -19*2  db/dec 


cp^  =  -11.75  deg/dec 


C’>>V 


^  J'i 

-  j  H 


V-V's 


Amplitude 
Ratio ,  db 


Therefore, 


'K 


jO.4^0 

19-2  .  11 .7^ 

20  ~  ^  132 


0.46?  ^  264.7  deg 


This  is  an  amplitude  error  of  5  percent  and  an  angular  error  of  less  than  1  deg. 

The  Siggy  of  Fig.  A-1  can  also  be  used  with  uhe  root  decoiBposition  technique 
to  estimate  A3,  due  to  a  gain  change.  For  a  6-db  increase  in  gain,  AK/K  =  1, 
root  decomposition  gave 


^1 

=  0.10 

^2 

=  -0.05  -  jO.345 

Therefore  from  Eq  154 

q1 

=  0.10 

which  is  in  error  by  5  percent,  and 

Sg  4  .0.05  -  jO.345 
4  0.348  .jr  262  deg 

which  is  a  magnitude  error  of  29  percent  and  an  angle  error  of  lesc  than  5  (itg* 
4.  Method  Using  the  Closed-Loop  Bode  Asymptotes 


The  closed- loop  Bode  amplitude  asymptote  will  be  unity  (O  db)  for  frequencies 
up  to  0.637  and  at  that  point  it  will  break  down  at  -40  db/dec.  For  q2,  then. 


^asy  ^ 


-  4q  -  Xu  -  0 


and  from  Eq  1 62 


,a  ■  (0-637)^ 


«  =  2j(6.%0)  '  4  270  des 


This  is  an  error  of  8  percent  in  magnitude  and  5*5  deg  in  direction. 


From  Eq  1 60 


=  (5-10)  (0.01 56)  =  0.0796 


This  is  an  error  of  16  percent. 


B.  LOW  GAIN 


We  will  now  consider  the  case  when  the  gain  is  set  so  that  there  is  a  second- 
order  pole  on  the  negative  real  axis.  For  this  case  the  important  parameters  are 


^  (7V^  -  27)  = 
i  (6  +  2V^)  = 
1  (6  -  V^)  = 


=  1.128 


5.055 


0.472  (second-order  pole) 


I .  Direct  Calculation 


The  numerator  and  denominator  polynomials  and  their  derivatives 


=  s^  •«-  6s^  +  “^s 


=  +  12s  +  5 


=  6s  +  12 


98 


Then  using  Eq  1 30  for  the  niunerator  and  denominator  derivatives  method, 

sj  =  - — - -1  = 

Uc(3^a/Ss^)  +  _  -6^2  +  12 

s — q2 

=  -  ^  (^9  -  9V^)  =  -0.246 

Sj,  =  r - - 1  =  - _ 

Lk(^/Ss)  +  (Sp/^s)J  3q?  -  I2q,  +  5 

5 — q^  I  I 

=  ^  (7V^  -  27)  =  0.0557 

Note  that  for  this  case  the  evaliaation  of  the  gain  sensitivities  only  determine 
two  of  the  modal  response  coefficients, 

«i  =  4 

%  =  -4 

To  evaluate  the  third  coefficient,  ,  we  must  sake  additional  calculations. 

Eor  this  case  the  easiest  method  is  to  use  Eq  115,  i.e., 

Q,  +  02^  =0 

or  Qg^  =  -Q,  =  .si 

The  same  result  could  he  obtained  from  Eq  112,  which  for  this  case  reduces  to 


-%  _  -  ^  (^9  -  9V^) 

“  ^2  ^  (6  +  2V^)  -  i  (6  -  V^) 

=  -  m 


Fcr  the  second  method  of  direct  calculation,  called  the  summation  of  terms 
method,  from  Eq  157  we  have 


ik^) 


ikhf 


We  can  set  up  a  computation  procedure  similar  to  that  used  for  the  high-gain  case, 
i.e. , 


“m"  '«• 


Pk 

(Pk  ■  ^2^ 

k  ■  '=2) 

0 

-0.47247 

4.4797 

1 

0.52753 

3.595*^ 

5 

4.52753 

0.0488 

s:|i 

2  - 

— - - ^  =  8.1219 

'1.  - 

S‘t  =  -0.246 

K 


Pk  (Pk  -  ^1 ) 

(pk  -  ^1, 

0  -5*05505 

-0.198 

1  -4.05505 

-0.247 

5  -0.05505 

-18.165 

y(  ^  'i  =  -18.610 
^  VPk  - 

s)^  =  0.0557 

100 

■*' • 

•  ’*•"  *- 

2.  Root-Locus  Methods 


For  the  gain-perturbation  method  the  gain  can  be  increased  to  2.070,  for  which 
the  root  was  previously  determined  to  be  -0.450  +  j0.450.  Then  from  Eq  138 


S 


2 

K 


i.i28[o.450  -  jO.450  -  0.472]^ 

2.070  -  1 .128 

-0.242  +  jO.0237  =  0.243  .^174.4  deg 


Even  though  the  ratio  of  the  gains  is  more  than  1  .8,  the  estimate  is  within 
2  percent  of  the  correct  amplitude  and  6  deg  of  the  right  direction. 

Similarly  for  , 

1. 128[5. 099  -  5.055] 

sl  =  -  =  0.0526 

2.070  -  1  .128 

This  is  a  2-percent  error. 


Using  the  phase -perturbation  method  for  ^2^  ^  perturbation  along  a 

line  midway  between  branches  of  the  root  locus.  Selecting  a  perturb>-tion  of 
0.25  at  an  angle  of  45  deg  gives  a  phase  change  of  l4.6  deg.  Therefore,  from 


Eq  159 


-0.245 


This  is  an  error  of  less  than  1  percent. 

For  q^  a  perturbation  of  jO.1  was  used.  This  gives  a  phase  change  of  63-7  deg. 


Then, 


j57.5[-j0.i] 


0.090 


This  estimate  has  a  large  error,  68  percent,  but  this  should  not  be  surprising 
because  the  perturbation  selected  was  nearly  twice  as  large  as  the  distance  from 
q^  to  the  open-loop  pole  at  -5.  Under  these  conditions,  higher  order  effects  are 
certain  to  be  important. 


For  q2f  Eq.  15O  'becomes 


-2(2.50)  (20)  (0.2)^ 


s=-q2 


Two  methods  of  evaliiating  A2  were  used.  The  first  was  to  measure  the  slope  of 
the  Siggy  (this  is  A-j )  and  plot  this  slope  versus  frequency,  see  Fig.  A-5.  The 
slope  of  this  curve  is  then  A2*  The  value  measured  from  Fig.  A-5  was  75  dh/dec^. 
This  gives 

A  -0.281 


which  is  in  error  by  l4  percent. 

The  second  method  was  to  estimate  the  radius  of  curvature  of  the  Siggy  at  a 
frequency  of  0.472.  A  radius  of  50  db  appeared  to  match  the  curvature  of  Fig.  A-5 
The  vertical  scale  of  the  graph  is  20  db/in.,  ana  the  horizontal  scale  is 
0.4  dec/in.  Following  the  procedure  outlined  in  Section  IH  gives 

''s  =  (oTtS^)  30^  =  85.3»/dec2 

This  gives  a  sensitivity  of 


A  -0.247 

which  is  amazingly  close  to  the  exact  value  (less  than  1  percent  error).  The 
surprising  agreement  is  obviously  fortuitous  as  the  method  is  certainly  not  that 
accurate . 

For  q.]  we  are  again  faced  with  the  problem  of  being  on  an  extremely  steep 
portion  of  the  Siggy.  Rather  than  attempt  to  measure  the  slope,  we  will  estimate 
it  the  same  way  we  did  for  the  high-gain  case.  The  two  low  frequency  roots 
contribute  -40  db/dec  and  the  root  at  -5  contributes 

-20^5^1^)  =  -I838  dl>/dec 


105 


Therefore 


=  -1878  db/dec 


and 


ol  i  -20(^.055) 
ItBtS 


0.0539 


which  is  within  1  percent. 

For  the  gain-perturbation  method  a  gain  change  of  -6  db  {M./K  =  -1  /2)  was 
selected.  It  can  be  seen  from  Fig.  A-3  that  this  change  splits  the  second-order 
pole  into  two  first-order  poles  with 


Thus, 


q  =  0.1 3,  0.83 


!0.13  -  0.47  =  -0.3^ 

or 

0.83  -  0.47  =  0.56 


Either  value  of  Ag.2  could  be  used  in  Eq  154,  but  because  neither  value  has  a 
theoretical  advantage  the  average  magnitude  of  the  two  was  used.  Therefore, 


s2  i  =  .0.21,5 

■  2 

This  is  less  than  1  percent  error,  although  the  method  is  obviously  not  generally 
that  accurate. 

For  q^  we  make  use  of  the  fact  that 

1  k 

so  that  with  the  gain  change 


Then 


q^  =  6  -  (0.13  +  0.85)  =  5-04 


1  ^  5.04  -  5.055 
'k  ~ 

"  2 


0.05 


This  is  an  error  of  44  percent. 


4.  Method  Using  the  Closed-Loop  Bode  Asymptotes 


The  closed- loop  Bode  amplitiide  asymptote  will  be  imity  (O  db)  for  frequencies 
less  than  0.472  and  at  that  point  it  will  break  down  at  -40  db/dec.  For  q2,  then 


A  =  1 
asy 


“  ^q  ~ 


and  from  Eq  l60 


"/c 


■(0.472)^  =  -0.223 


This  estimate  is  within  10  percent  of  the  exact  value. 
For  q^ 


A 


asy 


/o.472\ 

V5.055/ 


=  0.00875 


X  =  x„  =  0,  X,  =  2 


and 


sj,  =  (5.055)  (0.00875)  =  0.044l 


which  is  an  error  of  18  percent. 
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